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ABSTRACT 


A  full-wave  solution  has  been  obtained  for  radio  wave  propagation  in  the 
presence  of  an  elevated  tropospheric  layer.  This  analysis  wTas  performed 
as  part  of  the  CNI  (Communication,  Navigation,  Identification)  experimental 
program.  The  tropospheric  layer  is  modeled  as  a  trilinear  refractive 
index  profile  with  a  sufficient  lapse  rate  so  as  to  result  in  an  elevated  duct. 
The  analytical  solutions  give  the  received  signal  level  for  air-to-air  propa¬ 
gation  paths  at  UHF,  the  corresponding  signal  fading  level,  and  the  space 
diversity  distance  to  insure  good  quality  reception.  The  analysis  performed 
here  can  be  applied  to  ionospheric  propagation  and  to  the  underwater 
acoustic  channel. 
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SECTION  1 
INTRODUCTION 

1. 1  Motivation  for  the  Problem 

In  the  course  of  designing  a  propagation  experiment  for  the 
Communication,  Navigation  and  Identification  (CNI)  Program  [1]  ,  for 

air-to-air  and  air-to-ground  communications  links  at  UHF ,  it  was  recog¬ 
nized  that  anomalous  tropospheric  conditions  could  lead  to  serious  de¬ 
gradation  (or  enhancement)  in  signal  strength  and  to  multipath  (inter¬ 
ference)  effects  similar  to  those  produced  by  the  earth.  These  facts 

prompted  an  initial  investigation  into  the  subject  of  tropospheric 

kk 

propagation.  It  was  learned  that  tropospheric  layers  in  the 

first  few  kilometers  above  ground  are  relatively  common,  but  their 
precise  effect  on  signal  strength,  fading  characteristics,  Doppler 
spectrum,  available  bandwidth,  space  diversity  and  relative  terminal 
heights  has  been  investigated  in  relatively  few  cases. 

Propagation  in  an  anomalous  atmosphere  has  been  studied  exten¬ 
sively  using  the  theory  of  geometrical  optics.  The  theory  was  quite 
successful  in  explaining  the  gross  features  of  the  results  of  air-to- 
air  propagation  tests  at  UHF  [2].  However,  the  theory  is  approximate 
and  it  suffers  from  the  drawback  of  all  geometrical  optics  treatments 
in  that  the  solution  breaks  down  at  the  caustics,  and  furthermore, 
cannot  predict  field  strengths  beyond  the  horizon.  Analytical  solu¬ 
tions  or  full  wave  solutions  for  elevated  layer  propagation  have  re¬ 
ceived  only  a  few  approximate  and  incomplete  treatments.  Furry  f5]  has 
obtained  approximate  solutions  which  are  valid  for  the  trapped  or  ducted 
waves  only,  Northover  [6]  has  treated  (at  VHF)  the  elevated  layer  as  a  dis- 

*  A  bibliography  is  listed  at  the  end  of  the  report. 

**  The  term  layer  means  a  stratum,  of  vertical  thickness  from  a  few 
meters  to  several  hundred  meters,  within  which  the  mean  vertical 
gradient  and/or  the  variance  of  refractive  index  are  much  greater 
than  elsewhere  (see  Section  1.2). 
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continuity  of  the  index  of  refraction,  which  is  an  unrealistic  model 
in  the  general  case,  and  Chang  [5 1  considered  (at  VHF)  only  the  trapped 
waves,  essentially  ignoring  the  effects  of  the  earth. 

As  a  result  of  this  initial  investigation,  it  was  seen  that  the 
present  insufficient  and  incomplete  knowledge  of  wave  propagation  at 
UHF  in  the  presence  of  elevated  layers  could  not  lead  to  the  effective 
planning  of  a  propagation  experiment.  Moreover,  only  more  precise 
analytical  results  along  with  a  more  realistic  model  of  the  layer  would 
lead  to  a  better  interpretation  of  the  final  test  data  from  such  an 
experiment.  It  is  with  this  objective  in  mind,  that  the  present  work 
was  undertaken. 

1. 2  A  Discussion  of  the  Problem 

We  will  give  in  this  section  an  outline  of  the  problem  and  the 
results  we  hope  to  obtain.  We  will  discuss  the  assumptions  that  we 
are  making,  and  delineate  the  part  of  the  work  that  is  original.  We 
will  also  include  a  description  of  the  troposphere  and  the  effects  of 
elevated  layers  on  received  signal  strength. 

The  problem  considered  is  that  of  determining  the  electromagnetic 
field  produced  by  a  vertical  magnetic  dipole  (transmitter)  in  a  radially 
inhomogeneous  medium  (troposphere)  on  or  above  a  spherical  earth.  A 
full  wave  solution  is  obtained  for  the  case  when  the  variation  of  the 
index  of  refraction  of  the  atmosphere  is  such  as  to  produce  an  elevated 
tropospheric  duct.  The  solutions  for  the  field  intensity  are  obtained 
as  a  function  of  frequency,  layer  height  and  vertical  thickness,  lapse 
rate  or  change  of  index  of  refraction  with  height  (dN/dh)  within  the 
layer,  and  as  function  of  transmitter  and  receiver  heights  and  separa¬ 
tions.  In  the  solution  of  the  problem  the  WKB  or  phase  integral  solu¬ 
tions  for  a  second  order  ordinary  differential  equation  are  used.  The 
work  of  previous  authors  is  extended  in  three  important  ways.  First,  we 
use  the  "extended"  WKB  solutions  of  Langer  [  6  ]  .  This  approach  allows 
trie  closing  of  rhe  "gap"  in  the  regular  WKB  solutions.  It  involves  Hankel 
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functions  of  order  one-third.  Second,  the  Langer  solutions  valid  only 
for  a  single  transition  point  are  extended  to  be  valid  for  three  tran¬ 
sition  points.  This  is  essential  to  obtain  results  valid  in  the  elevated 
duct  and  it  is  accomplished  by  obtaining  solutions  in  three  regions  and 
joining  these  at  the  bottom  and  top  of  the  elevated  layer  heights.  And 
third,  the  resultant  solutions  for  the  trapped  and  leaky  waves,  are  ob¬ 
tained  from  the  solution  of  a  single  transcendental  equation.  This  re¬ 
sult  is,  to  the  best  of  the  author's  knowledge,  entirely  original. 

The  refractive  index  profile  of  the  troposphere  is  given  by  three 

linear  line  segments  of  different  slope.  While  at  first  sight  this  choice 

might  seem  artificial,  it  has  been  shown  by  ref ractometer  measurements  that 

index  of  refraction  profiles  do  take  on  such  abrupt  changes  and,  hence, 

the  trilinear  profile  becomes  a  valid  approximation.  In  reality,  all  that 

is  required  for  the  presented  solutions  to  apply  is  that  the  profile  be 

★ 

"sufficiently"  linear  in  the  three  segments  .  The  second  assumption 

we  are  making  is  to  replace  the  spherical  earth  by  a  flat  earth  and  modify 

the  index  of  refraction,  so  called  earth  flattening  approximation.  In 

this  case  the  solutions  are  valid  to  within  two  percent  up  to  ranges  of 

about  half  the  radius  of  the  earth,  and  this  independent  of  the  frequency. 

On  the  other  hand,  the  error  in  the  height-gain  functions  increases  with 
5/2 

elevation  as  h  and  is  proportional  to  the  frequency  f  7  ]  .  At  a 
frequency  of  300  MHz  and  an  elevation  of  5000  feet,  for  example,  the 
height-gain  function  is  only  in  error  by  less  than  two  percent.  The 
third  assumption  is  that  the  earth  is  a  perfect  conductor.  This 
is  a  good  approximation  for  horizontally  polarized  waves  and  for  the  low 
elevation  angles  of  interest  here. 

We  will  now  turn  to  a  discussion  of  the  troposphere  and  the  effects 
of  elevated  ducts  on  communication  links. 


The  sufficiency  condition  is  given  in  Section  1.3. 
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The  elevated  inversion  layer  or  elevated  duct  is  characterized  by 
a  rapid  decrease  in  the  index  of  refraction,  n  over  a  relatively  short 
height  interval.  The  index  of  refraction  n  can  be  derived  from  basic 
thermodynamic  laws  and  can  be  shown  to  give 


n 


1  + 


77.6 


(P 


+  4810e ) 


where  p  is  the  atmospheric  pressure  (millibars) ,  T  is  the  temperature 
(absolute  scale)  and  e  is  the  partial  pressure  of  water  vapor  (milli¬ 
bars).  The  refractive  index  of  the  troposphere  does  not  to  a  first 
order  approximation  depend  on  the  frequency  for  wavelengths  longer  than 
1  cm.  For  waves  of  the  millimeter  range  considerable  losses  appear, 
which  can  be  considered  by  means  of  the  introduction  of  a  complex 
dielectric  constant  of  air.  Then  the  index  of  refraction  n  will  depend 
on  frequency.  In  this  report  we  are  not  interested  in  frequencies  higher 
than  10  GHz . 

In  view  of  the  small  departure  of  n  from  unity  (on  the  order  of 

—A  6 

10  )  it  is  often  convenient  to  use  the  quantity  N  =  (n  -  1)  10  ,  called 

ref ractivity .  Near  the  earth’s  surface  the  index  depends  on  the  climatic 

and  meteorological  conditions  and  varies  within  260-460  N-units.  On 

the  average,  the  value  of  N  changes  linearly  with  height,  and  for  the 

middle  latitudes  the  gradient  of  the  change  N  with  height  h  is 

HN 

—  ~  -40  N-units/km  . 

dn 

An  elevated  layer  is  shown  qualitatively  in  Figure  1.1  at  an  altitude 
between  1000  and  1100  meters. 
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FIGURE  1.1  ELEVATED  DUCT  — 
N-PROFILE 


Most  frequently  these  layers  are  caused  by  a  mass  of  warm  dry  air  over- 

lying  a  cool  moist  air  mass.  The  usual  occurrence  of  this  phenomenon 

is  at  heights  of  approximately  1  to  2  kilometers.  When  the  refractivity 

gradient  exceeds  157  N-units/km,  within  the  layer,  the  layer 

becomes  an  elevated  duct.  The  changes  in  N  have  been  measured  to  be  as 
3 

high  as  5  x  10  N-units  per  km.  Sometimes  20-25  N-unit  changes  have 
been  measured  across  a  layer  having  a  thickness  of  a  few  meters  only. 

Even  though,  the  modified  index  or  earth  flattening  procedure  has 
certain  limitations  as  pointed  out  above,  it  offers  great  simplification 
and  a  large  insight  into  duct  propagation.  For  this  reason,  we  will  use 
the  concept  of  the  modified  index  of  refraction  to  give  a  qualitative 
description  of  propagation  in  elevated  ducts.  The  modified  index  of 
refraction  seeks  to  accomplish  two  aims:  (1)  its  vertical  gradient 
should  vanish  at  any  elevation  for  which  the  path  of  a  ray  launched 
horizontally  is  a  circular  arc  concentric  with  the  surface  of  the  earth; 


5 


and  (2)  it  should  be  readily  calculable  from  the  refractivity  N. 
For  this  reason  the  modified  index  is  defined  by  [8], 


M(h)  =  N  +  106  - 
a 


where  h  is  the  height  above  the  earth  and  a  is  the  radius  of  the  earth. 


FIGURE  1.2  ELEVATED  DUCT  — 
M  -  PROFILE 
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In  a  distribution  in  which  M  is  constant  rays  are  straight,  when  dM/dh 
is  negative  rays  curve  downward  and  when  dM/dh  is  positive,  the  rays 
are  refracted  upward.  The  net  result  is  that  through  the  use  of  M(h) 
a  spherical  coordinate  problem  is  more  simply  analyzed  in  rectangular 
coordinates,  where  plane  wave  modes  have  replaced  spherical-harmonic 
modes.  An  M-profile  roughly  corresponding  to  the  N-profile  of  Figure 
1.1  is  illustrated.  We  will  use  a  profile  similar  to  Figure  1.2  in  the 
solution  of  our  problem. 

The  disturbing  effect  of  tropospheric  ducts  on  air-to-air  and  air- 
to-ground  communication  links  has  been  investigated  in  several  important 
experiments  [9,  10,  111.  In  general  three  types  of  disturbances 

have  been  observed  which  cannot  be  explained  in  terms  of  ground 
reflections  and  which  are  evidently  associated  with  the  tropospheric 
propagation  medium.  These  disturbances  are  (i)  intervals  of  range  in 
which  average  signal  has  decreased  15  or  20  dB  (radio-holes) ,  (ii) 
intervals  of  range  in  which  the  average  signal  exceeds  the  free-space 
value  by  10  or  12  dB  (anti-holes)  and  (iii)  intervals  in  range  in 
which  interference  type  of  lobes  (multipath)  of  large  amplitude  (10- 
20  dB)  appear.  It  has  been  observed  that  radio  holes  and  anti-holes 
occurring  far  from  within  to  far  beyond  the  standard  horizon  are  very 
common  air-to-air.  For  instance,  all  of  the  50  air-to-air  propaga¬ 
tion  flights  made  by  Wright  Air  Development  Center  in  1951-1953  en¬ 
countered  radio  holes;  these  flights  were  made  under  fair  weather 
during  May-October  in  a  three  year  period  in  Ohio  [9].  Radio  holes 
far  within  the  standard  horizon  are  not  so  common  ground-to-air. 
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1* 3  Outline  of  Approach 


The  procedure  is  to  reduce  the  problem  of  determining  the  electric 
field  strength  of  a  dipole  anywhere  in  space  to  the  solution  of  an  in¬ 
homogeneous  Hertzian  potential  wave  equation  with  radial  and  angular 
coordinate  dependence.  The  solution  is  subject  to  the  boundary  condition 
at  the  earth,  the  Sommerfeld  radiation  condition,  the  continuity  of  the 
solution  plus  its  derivative  at  all  non-singular  points  and  the  condition 
that  the  solution  must  remain  finite  at  the  dipole  source.  By  separation 
of  variables  the  Hertzian  potential  is  represented  as  an  infinite  sum  of 
radial  and  angular  dependent  functions.  The  angular  functions  are 
Legendre  functions  and  the  radial  functions  are  as  yet  unspecified.  The 
coefficients  of  the  summation  can  be  obtained  from  the  boundary  conditions. 
The  resulting  series  is,  however,  a  very  slowly  converging  series.  To 
improve  the  rate  of  convergence  of  this  series  the  Watson  transformation 
is  applied.  The  individual  terms  of  this  resulting  series  are  evaluated 
at  the  roots  for  which  the  Wronskian  vanishes.  This  will  give  a  general 
solution. 

The  refractive  index  profile  is  then  chosen  as  to  give  three  trans¬ 
ition  points;  this  corresponds  to  an  elevated  duct.  The  radial  dependent 
functions  mentioned  above  are  evaluated  by  an  approximate  solution  to  the 
radial  wave  equation  with  the  "extended  WKB"  method  due  to  Langer  [6], 
at  each  transition  point  and  in  three  separate  regions.  The  solutions 
are  Bessel  functions.  The  solutions  in  the  three  regions  are  then  joined 
by  the  continuity  conditions.  The  fact  is  then  emphasized  that  the 
propagation  problem  reduces  to  the  solution  of  an  eigenvalue  problem  for 
an  ordinary  differential  equation.  Hence,  the  solution  has  to  satisfy 
the  boundary  condition  at  the  surface  of  the  earth  and  the  radiation 
condition  at  infinity.  The  resulting  set  of  complex  eigenvalues  and  the 
associated  eigen-functions  represent  a  complete  set  of  solutions  to  the 
problem.  Substituting  this  set  of  eigenvalues  and  eigen-functions  in 
the  general  series  solution,  we  can  determine  the  field  strength  due  to 
the  radiating  dipole  at  any  point  in  space. 


8 


In  the  solution  of  this  problem  it  will  be  seen  that  an  eigen¬ 
value  ordinary  differential  of  the  form  occurs: 

,2 

d  u  2 

— +  k  [y  (h )  +  A  ]  u  =0  (1.1) 

.,2  o  n  n 

an 

where  A^  and  u^  are  the  eigenvalues  and  eigenfunctions,  respectively, 

k  =  2tt / A ,  and  y(h)  is  related  to  the  modified  index  of  refraction  by 

y (h)  =  2  x  10  (M(h)  “  M  ).  The  functions  u  must  satisfy  the  boun- 

o  n 

dary  condition  at  the  earth  and  the  radiation  condition  at  infinity. 

The  form  of  the  solution  to  this  differential  equation  will  depend  upon 

the  term  in  brackets  and  in  particular  upon  the  eigenvalues  A^  which 

are  in  general  complex.  When  the  A^  have  imaginary  parts  (positive) 

which  are  comparable  in  magnitude  to  their  real  parts,  the  corresponding 

wave  is  said  to  be  "leaky".  This  is  because  the  value  of  the  function 

u^(h)  is  found  to  increase  rapidly  with  height.  When  the  imaginary 

part  of  A^  becomes  small  and  the  real  part  of  A^  is  comparable  to 

(-M(h^))  (see  Figure  1.1),  we  have  "transitional”  waves.  And  when 

the  imaginary  part  of  A^  ^  -M(h^),  the  waves  are  ’’trapped",  because  the 

functions  u  (h)  has  its  largest  value  between  h  and  h  .  We  will  use 
n  a  c 

the  "extended"  WKB  (Wentzel,  Kramers,  Brillouin)  solutions  of  equation 
(1.1)  in  the  three  cases  of  leaky,  transitional,  and  trapped  waves. 

In  general,  the  term  [y(h)  +  A  ]  vanishes  at  three  points  which, 
though  complex,  are  extremely  close  to  the  real  values  h^ ,  h^ ,  h^  shown 
in  Figure  1.1. 

Such  points  are  called  turning  points  or  transition  points.  This 
description  has  its  origin  in  quantum  mechanics  where  an  equation  very 
similar  to  equ.  (1.1),  namely  Schrodinger 1 s  equation,  has  been  the  sub¬ 
ject  of  considerable  analysis.  The  solutions  of  this  equation  make  a 
transition  from,  say  oscillatory  to  exponential  at  the  transition  point. 

The  Jeffreys,  Wentzel,  Kramers,  Brillouin  (JWKB)  or  WKB  method  is 
a  procedure  for  finding  approximate  solutions  to  a  second  order  differ- 
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M 


FIGURE  1.3  M-PROFILE 


ential  equation  of  the  form 


u"  +  f(z)u  =  0 


(1.2) 


in  some  region,  R,  in  which  f(z)  is  analytic  and  contains  no  zeros, 
and  in  which  f(z)  varies  sufficiently  slowly.  This  means  that  the 
first  two  derivatives  of  f  have  to  be  sufficiently  small  throughout 
the  region  in  question,  with  f  ^  0  therein  [4].  The  WKB  solution  of 

(1.2)  is  given  by 
1 


(1.3) 


where 


z 


(1.4) 


o 
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•k 

and  A  and  B  are  arbitrary  constants.  However,  the  WKB  approximation 
always  fails  sufficiently  near  a  zero  of  f(z).  A  form  of  the  solution 
of  (1.2)  which  is  valid  up  to  the  zero  of  f (z)  itself  was  given  by 
R.  E.  Langer  [6].  This  solution  has  the  form 


u 


A  H 


(1) 

1/3 


(F)  +  B  H 


(2) 

1/3 


(F) 


(1.5) 


where  F  =  F(z)  is  as  defined  above  and  A  and  B  are  arbitrary  constants, 
and  the  H^/^  are  Hankel  functions  of  the  first  and  second  kind.  This 
approximation  satisfies  the  equation 


ai  +2  ad  _i  c' 

f2  4  q2  2  Q  , 


=  0 


(1.6) 


where  Q  =  f 


Hence,  if  we  can  show  that 


2  sii  . 

4  Q2 


5 .  £  .  i£<( 

36  r2  2  Q 


(1.7) 


the  solution  in  (1.5)  is  a  very  good  approximation  of  the  original 
differential  equation  (1.2). 

A  limitation  of  the  Langer  WKB  solution  extension  is  that  the 
method  applies  only  to  a  single  transition  point.  Hence,  in  our 
approach  to  the  problem  we  shall  obtain  "extended"  WKB  solutions  of 
the  differential  equation  in  a  region  including  only  one  transition 
point.  These  solutions  are  then  "joined"  at  the  border  of  each  region 
by  employing  the  continuity  of  the  functions  and  their  derivatives , 
so  that  the  general  solution  is  valid  in  the  semi-infinite  interval 
extending  from  the  surface  of  the  earth  to  infinity. 

_ 

Northover  [5  ],  p.  37ff,  shows  that  for  the  WKB  solution  to  hold 
we  must  have  |f,/f3/2|  <<  1  and  f"/fl/2  <<  i. 
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1.4  Contents  of  the  Report 


We  will  give  here  a  brief  outline  of  the  remainder  of  this  report. 

In  Section  2,  we  formulate  the  problem  in  a  precise  mathematical  way. 

The  three  dimensional  vector  solutions  of  Maxwell fs  equations  are  re¬ 
duced  to  a  one  dimensional  vector  by  the  introduction  of  the  Hertz 
potential.  By  placing  the  transmitting  source  on  the  polar  axis,  ad¬ 
vantage  is  taken  of  the  symmetry  of  the  problem.  As  a  result,  the 
Hertz  potential  has  only  an  r  and  0  dependence.  The  electric  and 
magnetic  fields  can  then  be  written  in  terms  of  this  Hertz  potential. 

In  Section  3,  we  obtain  the  slowly  converging  series  solutions  for  the 
Hertz  potential.  By  separation  of  variables  this  series  can  be  written 
as  the  product  of  the  radial  dependent ,  r  and  angular  dependent ,  0 
functions.  The  angular  functions  are  Legendre  polynomials,  and  the 
radial  functions  are  as  yet  unspecified;  since  they  depend  on  the 
particular  index  of  refraction  profile  assumed.  In  Section  4,  we  apply 
the  Watson  transformation  to  obtain  a  more  rapidly  converging  solution. 

This  is  done  by  expressing  the  series  solution  for  the  potential  as  a 
contour  integral  in  the  complex  plane  enclosing  an  infinite  number  of 
singularities.  By  Cauchy’s  residue  theorem,  the  more  rapidly  converging 
series  for  the  potential  is  now  a  sum  over  the  residues  enclosed  by 
the  contour.  In  Section  5,  we  use  Langer’s  ’’extended"  WKB  solutions 
to  obtain  approximate  solutions  of  the  radial  dependent  functions.  And 
in  Section  6,  we  obtain  the  radial  dependent  solutions  for  three  trans¬ 
ition  points.  In  Section  7,  we  transform  these  radially  dependent  sol¬ 
utions  from  that  of  a  spherical  geometry  into  a  "flat  earth"  geometry. 

In  Section  8,  we  present  a  model  for  the  layered  troposphere.  And  in 
Section  9,  we  obtain  the  expression  for  the  Hertz  potential  and  the 
electric  field.  In  Section  10  we  give  the  numerical  results  for  the 
potential.  The  discussion  in  this  section  includes  details  of  the  numerical 
evaluations.  Section  11  contains  a  summary.  Relegated  to  the  Appendices 
are  the  mathematical  details  of  the  contour  integral  representation  and 
the  modified  Hankel  functions.  Also  included  there  is  the  justification 
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for  representing  the  three  dimensional  vector  solutions  of  Maxwells 
equations  by  the  one  dimensional  Hertz  potential. 
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SECTION  2 


PRECISE  MATHEMATICAL  FORMULATION  OF  THE  PROBLEM* 


We  assume  the  earth  to  be  a  sphere  of  radius  a,  and  introduce  spher¬ 
ical  polar  coordinates  r,  0,<P  as  defined  in  Figure  2.1.  We  are  given 
a  vertical  magnetic  dipole  of  moment  ml  with  a  harmonic  time  var¬ 
iation  e  at  a  point  D(b,  0,  0)  and  we  are  required  to  determine 

the  electromagnetic  field  at  an  arbitrary  point  P(r,  0,  9)  on  or  above 

2  2  1/2 

the  surface  of  the  earth  and  at  a  distance  R  =  (r  4*  b  -  2br  cos  0) 

from  D.  We  assume  the  earth  to  be  immersed  in  an  inhomogeneous  medium 

whose  dielectric  constant  k  is  a  function  of  radial  distance  r  only, 

k  =  k(r);  all  effects  of  the  ionosphere  have  been  neglected  .  Inside 

the  earth  k(r)  is  assumed  to  be  a  constant,  k^.  In  any  medium  the 

propagation  constant  is  i<  =  [irny  (a  4*  ime)]!/^  so  that  in  air  k  =  2tt/X, 

1/2  ° 

while  in  the  earth  k9  =  k  (K  -  60  icrX)  ,  where  K  is  the  relative 

*  o 

dielectic  constant,  a  is  the  conductivity  (in  mhos  per  meter)  and  the 

permeability  y  is  assumed  to  be  the  same  inside  and  outside  the  earth. 

2  2 

In  this  problem  k  (r)  =  m  ye(r)  for  the  troposphere. 

The  electromagnetic  field  must  satisfy  Maxwell’s  equations,  both  in¬ 
side  and  outside  the  earth,  it  must  possess  at  the  point  D  a  singularity 
corresponding  to  that  of  the  dipole,  the  electric  field  and  magnetic 
field  must  be  continuous  everywhere,  while  at  infinity  the  field  vectors 
must  satisfy  the  Soramerfeld  radiation  condition. 

Let  E  and  H  represent  the  electric  and  magnetic  field  vectors, 
respectively;  the  factor  exp  (imt)  being  understood  throughout.  Then 
Maxwell’s  equations  in  m.k.s  units,  become  the  following: 


The  material  in  this  and  the  following  sections  is  part  of  the  author’s 
dissertation  for  the  Ph.D.  at  the  University  of  Pennsylvania, 
Philadelphia. 

**  It  has  been  shown  experimentally  that  at  frequencies  higher  than 
30  MHz  the  ionosphere  has  negligible  effects  on  terrestrial  wave 
propagation. 
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Figure  2.1  PROBLEM  GEOMETRY 


V  x  E  =  -iaiyH 


V  •  H  =  0 


(2.1) 


V  x  H  =  iojeE 


V  •  eE  =  0. 


It  has  been  shown  by  Friedman  [ 12  ]  that  in  spherical  coordinates  the 
electromagnetic  field  given  by  (2.1)  can  be  obtained  from  a  Hertz 
vector  II  which  has  only  its  radial  component  non-zero,  even  when  the 
dielectric  constant  is  a  function  of  r.  Since  the  dipole  is  located 
on  the  polar  axis,  symmetry  exists,  and  the  vector  II  will  not  depend 
on  V  but  only  on  the  coordinates  r,0.  We  define  II  =  (rU(r,0)  ,0,0) . 
In  terms  of  U  the  electric  and  magnetic  fields  are: 


V0’  V0’  v° 

H  =  -  1  32(rU) 

0  r  3r  30 


_  ,  3U 

E<0  -  30 


(2.2) 


H  =  -(k2  +  — j)  (rU). 
3r 


Equations  (2.2)  represent  the  field  of  a  magnetic  dipole  and  are  sat¬ 
isfied  for 


U 


-ikR 


M 


e _ 

4tvR 


(2.3) 


where  M  is  proportional  to  the  magnetic  moment  of  the  dipole,  its 
value  being 

M  =  (2.4) 

where  I  is  the  electric  current  in  the  loop  of  area  S  and  b  is  the 
distance  from  the  center  of  the  earth  to  the  dipole. 
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The  scalar  wave  function  U(r,0)  roust  satisfy  several  conditions. 
First  U(r,0)  must  satisfy  the  partial  differential  equation 

I.  V2U  +  k2u  =0  (2.5) 


in  the  exterior  region  r  >  a,  except  at  the  dipole,  and  in  the  interior 
region  r  <  a.  Second, 


II.  u  and  ■—  (rU) 
3r 


(2.6) 


roust  be  continuous  everywhere  except  at  the  source.  Thirdly,  as  R  ->  0 


III. 


U  M 


-ik  (b)R 
e 

4ttR 


(2.7) 


This  means  that  in  the  neighborhood  of  the  dipole  the  potential  U  be¬ 
comes  that  of  an  oscillating  dipole  in  an  homogeneous  medium  whose 
permittivity  is  equal  to  the  local  value  at  the  dipole  position. 

Finally,  U  must  satisfy  the  Sommerfeld  radiation  condition  given  by 

IV.  lim  r  (~  +  ik  U)  =  0  .  (2.8) 

3r  o 

^•-x» 

The  fact  that  U(r,0)  possesses  a  singularity  at  the  dipole,  r  =  b, 

0=0,  while  otherwise  it  satisfies  the  wave  equation  (2.5),  can  be 
expressed  very  concisely  using  the  Dirac  6-function.  We  write 


V2U  +  k2„  ,  _M  6  (r-b)6 (8) 
2Trr  sin  0 


(2.9) 


where  M  is  the  source  strength  and  the  factor  27Tr  sin  0  is  designed 
to  normalize  the  solution. 
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SECTION  3 


SERIES  SOLUTION  OF  THE  SPHERICAL  HERTZIAN  WAVE  EQUATION 


We  shall  now  obtain  a  general  solution  to  the  inhomogeneous 
Hertzian  potential  wave  equation  given  by  (2.9),  subject  to  the  re¬ 
quired  boundary  conditions  discussed  in  the  previous  section.  Rather 
than  writing  the  total  field  as  the  sum  of  a  primary  field  due  to 
the  dipole  source  and  a  scattered  field  due  to  the  earth,  we  shall  use 
a  self-consistent  method.  In  particular,  we  will  show  that  the  solu¬ 
tion  to  the  inhomogeneous  wave  equation  (2.9)  has  the  form  of  an  ex¬ 
pansion  in  zonal  harmonics, 


U(r,0) 


u  (r)  P 
n  n 


(cos  0) 


(3.1) 


where  n  are  integer  values 


and  u  (r) 
n 


are  solutions  to  the  equation 


d  (ru  ) 


dr 


*U[, 


k  (r)  - 


n(n+l) 


(ru  )  =  6  (r-b), 

v  n  2irr  v  * 


(3.2) 


or 


d  un  ,2  dun  .  f,  2,  „  n(n+l)l  -M  „  ,  ,  , 

"^2  +  7  “d?  +  [k  (r)  "  -7— J  Un  =  2^7  6  (r"b) 


(3.3) 


To  show  this  we  use  the  orthogonality  property  of  the  Legendre  poly¬ 
nomials  P  (cos  0)  in  (3.1)  and  find  that 
n 

r 

u  (r)  =  /  U(r,e)  Pn  (cos  0)  sin  0  d0  .  (3.4) 

0 

If  we  multiply  equation  (2.9)  (for  spherical  coordinate  system)  by 
Pn(cos  0)  sin  0  d0  and  integrate  between  0  and  it,  this  equation  transforms 
into  the  ordinary  differential  equation  (3.2). 
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The  foregoing  implies  that  for  U(r,0)  given  by  (3.1)  to  be  a 
solution  of  the  inhomogeneous  wave  equation  (2.9),  the  coefficients 
u^(r)  must  satisfy  eq.  (3.2).  Furthermore,  since  the  Legendre  poly¬ 
nomials  form  a  complete  set,  U(r,8)  as  given  by  (3.1)  is  uniquely 
specified.  Hence,  U(r,0)  given  by  eq .  (3.1)  is  the  total  solution 
of  the  inhomogeneous  wave  equation. 

The  boundary  conditions  that  had  been  imposed  on  U  in  Section  2 
must  now  be  satisfied  by  (ruR) ,  namely 

1.  (ru  )  and  (ru  )  must  be  continuous 

n  dr  n 

2.  u^  must  be  regular  at  r  =  0 

3.  The  Sommerfeld  radiation  condition 

lim  1 (ru  )  -  ik  (ru  )  '  =  0 
1  dr  n  on1 

must  be  satisfied. 

In  order  to  solve  equation  (3.2)  we  must  know  the  solutions  of 
the  homogeneous  differential  equation 


dr 


2  (tun>  + 


k2  (r)  - 

r 


(ru  )  =  0 
n 


(3.3) 


When  k  (r)  =  k^  ,  a  constant  (case  of  airless  atmosphere),  the  solutions 

of  (3.5)  are  linear  combinations  of  the  Hankel  functions  of  half-integer 

(1)  (2)  2 

order,  namely  Hn4.i/2^r^  anc*  n-fl/2^*  ^  ^  (r)  is  not  a  constant>  ttie 

solutions  of  (3.5)  behave  like  Bessel  functions,  but  their  exact  behavior 
will  depend  on  the  variation  of  k(r)  with  r. 

3. 1  General  Solution  for  the  Radial  Functions 

We  let 

u  (r)  =  a  V  (r)  for  r  <  b 

n  n  n  — 

u  (r)  =  b  Q  (r)  for  r  >  b 

n  nn  — 


19 


where  the  V^(r)  and  Q^(r)  are  independent  solutions  of  the  homogeneous 
differential  equation  in  (3.5).  They,  furthermore,  satisfy  the  re¬ 
quired  boundary  condition  in  their  respective  domains.  The  general 
solution  (3.1)  is  then  given  by 


U(r,e)  =  P  (cos  0) 

^0  2  n 


a  V  (r) 
n  n 


b  Q  (r) 
n  xn 


r  <  b 

(3.6) 

r  >  b  . 


The  coefficients  a^  and  b^  may  be  determined  by  applying  the  required 
boundary  condition  at  r  =  b.  At  r  =  b,  the  functions  (ru^)  are  con¬ 
tinuous  and  their  derivatives  have  a  junp  of  magnitude  -M/2tt  b.  From 
the  definition  of  V^(r)  and  Q^(r)  above,  the  functions  u^(r)  are  re¬ 
gular  at  r  =  0  and  satisfies  the  radiation  condition.  Then  we  can 
prove  that 


a 

n 


-M 

2ttW 


n 


Qn(b) 


(3.7) 


b 

n 


-M 

2ttW 

n 


Vb) 


(3.8) 


where  W  is  the  Wronskian  of  (rQ  )  and  (rV  )  evaluated  at  r  =  b.  Sub- 
n  xn  n 

stituting  (3.7)  and  (3.8)  in  (3.6),  we  obtain 


U(r,0) 


oo 

zfi  E 

4ir  " 


n=0 


(2n+l) 

W 

n 


P  (cos  0) 


n 


Vb)  Vr> 

Qn(r)  Vn(b) 


r  1  b 

(3.9) 

r  >  b  . 


This  is  the  solution  to  our  problem  in  its  most  general  form  for  an 
arbitrary  variation  of  the  index  of  refraction  of  the  atmosphere. 
Equation  (3.9)  represents  Greens'  function  for  the  inhomogeneous  wave 
equation  (2.9).  The  general  requirement  of  the  reciprocity  of  Greens' 
function  is  satisfied  owing  to  the  fact  that  eq.  (3.9)  is  symmetric  in 
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r  and  b.  The  reciprocity  with  respect  to  angles  0  and  0^  (we  con¬ 
sidered  the  case  in  which  the  latter  is  zero)  can  be  expressed  by  re¬ 
placing  cos  0  by  cos  (0  -  0^)  which  is  symmetric  in  0  and  0q.  It  has 
been  shown  that  the  solution  (3.9)  can  be  used  for  numerical  computa¬ 
tion  only  if  ka  is  small  compared  to  unity.  For  example,  at  100  MHz 
about  a  million  terms  are  required  to  evaluate  the  series.  In  the 
following  sections  the  solution  (3.9)  will  be  transformed  so  that  it 
can  be  readily  computed  for  large  values  of  ka . 
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SECTION  4 


THE  WATSON  TRANSFORMATION 


The  series  expression  for  U(r,0)  is  very  slowly  converging.  For 
example,  for  ka  of  the  order  of  10^,  where  k  =  2tt/X  and  a  is  the  radius  of  the 
earth,  approximately  10^  terms  of  the  series  are  required  to  evaluate 
the  value  of  U(r,0)  at  one  point.  Such  numerical  calculations  are 
even  prohibitive  in  this  age  of  high  speed  computers.  A  way  around 
this  laborious  task  was  found  by  Watson  [13].  Watson  showed  that  the 
series  (3.9)  can  be  expressed  as  a  contour  integral  in  the  complex  plane 
enclosing  an  infinite  number  of  singularities.  By  Cauchy1 s  residue 
theorem  this  integral  is  then  equal  to  the  sum  of  the  residues  enclosed 
by  the  contour.  The  net  effect  of  the  Watson  transformation  is  to  have 
a  much  faster  converging  series.  The  difficulty  that  is  left,  is  to 
locate  the  singularities,  namely  zeros  of  the  Wronskian  (see  Appendix  B) . 

From  eq.  (B.12)  we  have  the  series  representation 


U(r,0) 


M 

2 


E 


v  +  —  (cos  (it  -  e)) 
sin  vir 


dW 
ds 


s=v 


Vb)  Vr) 

Qv(t )  Vv(b) 


r  _<  b 

r  b 
(4.1) 


where  v  represents  the  zeros  of  the  Wronskian  Wg  and  the  derivative 
(d/ds)  is  with  respect  to  the  complex  variable  s. 

The  residue 


v  +  J  Pv(cos  (tt  —  0 ) ) 

sin  vtt  _ _  9 

W 


Q  (b)  V  (r) 
v  v 


where 


s=v 


(4.2) 


(4.3) 
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can  be  interpreted  in  another  way.  Since  the  Wronskian  Wg  =  0 ,  the 
two  functions  Q  (r)  and  V^(r)  are  linearly  dependent,  that  is  V^(r)  = 
a  Q(r).  This  means  that  the  function  Q^(r),  besides  satisfying  the 
Sommerfeld  radiation  condition,  also  satisfies  the  condition  of  regu¬ 
larity  at  r  =  0.  It  follows  that  v  is  an  eigenvalue  and  Qv(r)  is  an  eigen¬ 
function  of  the  following  problem: 

Find  the  values  of  v  for  which  V(r),  the  solution  of  the  differ¬ 
ential  equation: 


d2(rV) 

dr2 


+ 


j  —1  (rV)  =  0 
r 


(4.4) 


is  regular  at  r  =  0  and  satisfies  the  radiation  condition  at  infinity, 
that  is 


lira  \^~  (rV)  -  ik  (rV) I  =  0  .  (4.5) 

dr  o 

17 -K*3 


It  will  be  seen  in  later  sections  that  the  zeros  of  the  Wronskian, 
v,  will  be  very  large,  of  the  order  of  magnitude  of  ka.  This  allows  us 
to  use  asymptotic  formulas  for  the  Legendre  functions  P  (cos  (it  -  8)). 

We  have 


P  (cos 

v 


e) 


0 


tt(v+1)  sin  0 


1/2 

) 


cos 


[  (V  +  j) 


(4.6) 


for  large  v,  so  that 

P  (-cos  8)  P  (cos  (ir  -  8)) 

V _  _ 

sin  vir  sin  irv 


-  ( 


2i 


1/2 


7T  (v  +  1)  sin  0 


) 


exp 


{-i 


(V  +  j)  e } 


(4.7) 


23 


Next,  using  the  fact  that  the  dielectric  constant  of  the  earth 
is  constant,  it  has  been  shown  [12]  that  the  boundary  condition  at 
the  surface  of  the  earth  can  be  written  as 


hsf  <rV  - r  <rvl  * 0 

L  J  r=a 


r 


1/2 


(4.8) 


where  and  k^  are  the  wave-numbers  of  the  earth  and  free  space,  res¬ 
pectively.  A  perfectly  conducting  earth  implies  =  00  and  hence 
(4.8)  reduces  to 


CrVs) 


r=a 


0  . 


Equation  (4.8)  can  also  be  written  in  a  preferred  form: 


(4.10) 


(rv 

V  (rvs)  , 


=  r 


(4.11) 


r=a 


We  shall  use  this  result  presently. 

It  is  shown  in  Appendix  C,  that  with  the  aid  of  equations  (4.10) 
and  (4.11),  the  expression  (4.1)  can  be  written  for  an  inhomogeneous 
earth  and  a  perfectly  conducting  earth,  respectively,  as 


u(r  0,  m  y  a*-.-!)  vb>  ! 

'  4a2  V  sin  ”  «v(a>  <*v(a)  r 


(cos  (it  -  0)) 


8t  rQ 


(4.12) 


r=a 

t=v 
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U(r,0) 


-M  V"  (2v+l) 

4br  sin  vtt 

v 


(bQv(b))(rQv(r)) 

37  "V] 


P  (cos  (tt  -  6)  ) 

b  (rQb)L 


a 
t=v 


(4.13) 


Using  the  approximation  (4.7),  we  obtain 


U(r,6) 


yr(2v+l)  il 
2  sin  0  J 

a  v 


1/2 


exp 


[i(v  +  p  e] 


and 


Vb)  Qv(r) 

Qv(a)  Qv(a) 


at 


rQ, 


J  r=a 
t=v 


(A. 14) 


U(r ,6)  =  -M  ^2 


(2v+l)  i 

tt  sin  0 


1/2 


exp  [i(v 


e] 


Vb>  Qv(r> 


37  <rV]  [at  Wt1] 

J  r=a  L  J  r= 


a 
t=v 


(4.15) 


Equation  (4.14)  and  (4.15)  are  the  expressions  for  the  Hertz  potential 
satisfying  the  boundary  condition  (4.11)  and  (4.10),  respectively. 

The  formulation  in  (4.14)  (also  (4.15))  has  the  advantage  that 
the  same  formula  for  U(r,0)  represents  r  >  b  and  r  <  b,  in  contrast 
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to  eq.  (4.1).  The  expression,  however,  breaks  down  for  6  =  0  since 
for  small  0, 


^  „  \  sin  vtt  , 

P  (-cos  0 )  ->  -  log  0 

whereas  the  function  P^  represents  is  to  be  regular  for  0=0  and 
r  ^  b.  It  implies  that  the  whole  "ray"  0=0  must  be  considered  a 
singularity  of  the  representation  in  (4.14)  and  (4.15).  This  condition, 
however,  presents  no  problem  in  this  investigation,  since  we  will 
always  be  interested  in  fields  well  removed  from  the  0=0  line. 
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SECTION  5 


THE  LANGER  APPROXIMATE  FORMULAS 


In  Section  4  it  was  shown  that  the  problem  of  determining  the 
Hertz  vector  U(r,0)  and  then  the  electromagnetic  field  depends  upon 
the  solution  of  the  eigenvalue  problem  for  the  differential  equation 
(4.4).  In  this  section  we  will  discuss  the  "extended11  WKB  approxi¬ 
mations  due  to  Langer  to  obtain  solutions  of  this  differential  equa¬ 
tion.  The  extended  form  is  valid  in  the  neighborhood  of  the  transition 
point  as  well  as  at  the  transition  point  itself. 

Langer  has  demonstrated  that  the  solutions  of  the  so  called  re¬ 
lated  second  order  ordinary  differential  equation  for  a  single  trans- 

2 

ition  point,  i.e.,  a  zero  of  [$  (z)  -  0)(z)]  in  eq .  (5.1), 

g(z)  +  [(f)2  (z)  -  0)(z>]  g(z)  =  0  (5.1) 


is  given  by 


g(z)  =  nz)  s1/3  (c) 


j  =  1,  2  (5.2) 


in  which  the  symbols  stand  for  the  Hankel 

th 

1/3,  of  the  j  kind,  and 


functions  of  order 


nz)  =  S1/6(z)  <f>  1/2(z) 


■f 

J  Z1 


r,  (z)  =  I  <j>(z)  dz 

Z1 

T* 1  (z) 


u(z)  = 


nz) 


where  z^  is  a  zero  of  [ 4s  (z)  -  w(z)] 


(5.3) 

(5.4) 

(5.5) 
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Although  the  formula  for  ¥ (z)  seems  to  assign  a  singularity  to  this 
function  at  z^,  it  is  found  that  this  singularity  is  removable,  and 
the  ^(z^)  is  then  different  from  zero.  The  function  oj(z)  is  thus 
bounded  in  a  region  about  z^. 

2 

If  it  can  be  shown  that  oo(z)  is  very  small  compared  to  <j>  (z) 
over  the  considered  region  of  z,  then  (5.2)  may  be  considered  as  a 
solution  of  eq.  (5.1)  to  practical  accuracy.  In  that  case 

g  (r)  =  rQv(r)  , 

and 

+2(r)  -  ko2  [k2(r)  -  ^*4]  • 

o 

In  the  region  of  the  zero,  if 

U2(r)|  «  1  ,  (5.8) 

then  the  condition 


(5.6) 

(5.7) 


|  to  (z  )  |  <<  1 


(5.9) 


must  be  satisfied  for  (5.6)  to  be  an  approximate  solution  to  the 
wave  equation.  If  these  conditions  are  satisfied  in  the  interval  of 
interest,  for  which  there  is  only  a  single  transition  point  at  r  =  r^ 
then  the  general  solution  is  given  by 


rQ 

v 


k  "1/3  M  +  Bv  HU3 


(5.10) 


where 


and  4)(r)  is  given  by  (5.7).  Thus  eq .  (5.10)  represents  the  approxi¬ 
mate  solution  in  the  region  that  includes  one  transition  point.  For 
the  regions  of  the  second  and  third  transition  point,  solutions  of  a 
similar  form  apply,  each  valid  only  in  a  region  including  one  transi¬ 
tion  point. 

2 

It  has  been  shown  that  if  <J>  is  linearly  varying,  aj(z..)  =  0  and 
hence  eq .  (5.9)  is  satisfied  as  well  as  the  condition  j  cu  (z )  I  <<  <J>  (z) 
for  all  z. 


Furthermore,  Langer  has  shown  that  for  large  £  the  asymptotic 
forms  of  the  Hankel  functions  may  be  used.  Hence 


Hl/3  (t) 


=  c 


elL 

1  rl/2 


H«>  (C) 


=  c 


2  ^1/2 


(5.12) 


(5.13) 


where  and 
these  functions 
obtain 


rQ 

v 


are  constants  that  depend  on  the  normalization  of 
If  (5.12)  and  (5.13)  are  substituted  in  (5.10),  we 


(5.14) 


It  may  be  recognized  that  (5.14)  are  the  nordinary  WKB,?  solutions. 

2 

This  form  of  the  solution,  however,  breaks  down  at  a  zero  of  $  . 

2 

For  the  discussion  of  conditions  on  the  permittivity  when  <f> 
is  not  varying  linearly  the  reader  is  referred  to  Doviak  and 
Goldhirsh  ‘[14]  and  Bremmer  [15]. 
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SECTION  6 


SOLUTIONS  FOR  THREE  TRANSITION  POINTS 


As  we  have  mentioned  in  the  previous  section  the  extended  WKB 
solutions  are  valid  in  a  region  including  only  one  transition  point. 
It  was  pointed  out  in  the  Introduction  that  for  the  problem  of  the 
elevated  duct,  three  transition  points  exist.  Hence  we  must  obtain 
three  separate  solutions  for  each  of  the  three  regions.  These  solu¬ 
tions  are  then  matched  at  points  between  the  transition  points  using 
the  continuity  boundary  conditions. 


We  will  represent  the  solutions  (rQ^)  to  the  homogeneous  wave 
equation  in  regions  I,  II  and  III  by  (rQ^)^,  (rQ  and  (rQ^)^,  res¬ 
pectively.  Region  I  extends  from  r  =  a  to  r  =  region  II  for 

r\  <  r  <  n0,  and  region  III  from  r  =  ru  to  infinity.  Since  (rQ  )3 
i  -  -  z  z  v  (2) 

must  satisfy  the  Sommerfeld  radiation  condition,  only  solutions 
are  allowed.  Hence,  we  have  in  Region  I: 


1/2 


1/2 


<rVl  '  Av  HU3  <‘l>  +  Bv  <*7>  H{/3  «!>■•  ‘6'1) 


in  region  II: 


<rV2  ’  cv  V 


1/2  1/2 

HW  <‘2>  +  Dv  <^>  "l/l  <C2);  (6'2> 


in  region  III: 


1/2 

’3,  (2) 


<rV3  =  Fv  Hi/3  <*3>  5 


(6.3) 


where 


(r)  dr 
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(rQ  ) .  and 
v  J 

pect  to  r) 
A  ,  B  ,  C 

\>  V  V 


where 


! 


■f 


?2  (r)  =  /  <Mr)  dr 


(6.5) 


?3(r)  =/  <Kr>  dr 


(rQ  ) .  (j  =  1,2,3)  (prime  means  differentiation  with  res- 
must  be  continuous  at  r  =  and  r  =  r\  m  The  coct  f  icit-n ts 

and  D  may  be  expressed  only  in  terms  of  F^.  Further.  1  r 


1/2 

x(1)  - 

mi 

Hd) 

H 1  /  3 

(W 

(ft.  7) 

1/2 

x(2)  = 

lnl 

u  (2) 

Hl/3 

(hnl} 

(6.3) 

1/2 

X(1)  - 
2nl 

C'1) 

^21 

Hd) 

Hl/3 

^2nl^ 

( h  .  4  ) 

<l2nl> 

^21 

1/2 

x(2)  - 
A2ni 

h(2) 

Hl/3 

(?2nl) 

(ft. 10) 

1/2 

x(1)  = 

2n2 

<^> 

*22 

h(D 

H 1  /  3 

^2n2^ 

(ft. 11) 

1/2 

x(2)  = 

X2n2 

<'2"2> 

^22 

h(2) 

Hl/3 

(r’2n2} 

(ft.l?) 

1/2 

x(2)  = 

3n2 

C3"2) 

*32 

h(2) 

Hl/3 

^3n2^ 

(ft.  13) 

’lnl 


f 

=  I  4>dr 
«/r 


(6.1-0 
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(6.15) 


(6.16) 


(6.17) 


Using  the  above  in  equations  (6.1)  to  (6.3)  and  applying  continuity 
conditions  at  n^  and  we  have  the  set  of  four  simultaneous  equations 


A 

a}1} 

+  B  X<2>  = 

c  + 

D  A«> 

V 

lnl 

v  lnl 

v  2ql 

v  2nl 

A 

X(D, 

+  B  X<2>' 

=  C  A<\>« 

V 

lnl 

v  lnl 

v  2ql 

+  DvX 

c 

x™ 

+  D  X<2>  = 

F  X<2> 

V 

2n2 

v  2q2 

v  3q2 

C 

x<x>' 

+  D  X^2^ * 

=  F  X  ^2« ' 

V 

2n2 

v  2n2 

v  3n2 

(2), 


(6.18) 

(6.19) 

(6.20) 
(6.21) 


which  gives  for  A  and  B 

v  v 


A  =  K  F 
v  vi  v 


B  =  K  F 
v  v2  V 


(6.22) 

(6.23) 
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where 


+K 


r.  (2)  (2)  (2), 

lnl  3r,2  L1  lnl 


a<2>-  l2  +  a<2>’ 

3n2  2  lnl 


x(2)  1 

ALL- 


,(2)  (2), 

lnl  3n2  t 


D 


J 

(6.24) 


where 


D  = 


,(1)  ,(2),  _  ,(1),  y  (2)  ^(1)  ^  (2)  ,  _  >  (1)  ,  y  (2) 


lnl  Alnl  Alnl  Alnl  X2n2  *2n2  A2n2  '‘2n2j 


(6.25) 


and 


X(1)  X*2*  L  +  X  '  X  '  L  +  A  (1)'  X*2*  L  +  >!’]  A  ^ \ '  L 
lnl  3n2  1  lnl  3n2  _2 _ lnl  3n2  5  hi  Jn^  4 


Lb. 26) 


where 


L  = 

,  (i), 

X<2>’ 

-  A^’  *<”• 

(6. 

.27) 

1 

2nl 

2n2 

2nl  2n2 

L2  = 

x(1) 

2nl 

x (2)  _ 
X2n2 

X(1)  XU) 

2nl  2n2 

(6, 

.28) 

L3  " 

x(2) 

A2nl 

,  (1), 
A2n2 

,  (1)  y(2), 

A2nl  2n2 

(6. 

.29) 

L4  " 

,  (2), 
A2nl 

A1) 

2n2 

,  (1)  i  ,(2) 

2nl  2n2  * 

("•  ■ 

.  30) 

The  functions  L^, 

l2,  1 

and 

L,  depend  only  on  the  layer  between 

''  L 

and  n0.  It  will  be  seen  later  that  this  interval  <  r  ^  n.;  also 
corresponds  to  the  tropospheric  layer. 

Substituting  (6.24)  and  (6.26)  into  (6.1)  gives 
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(rQ  )  =  F 

v  v 


C  1/2  c  1/2 

Kv1  (?7>  “i/3  (tl>  +  Kv2  (^>  Hl/3  «1> 


(6.31) 

Substituting  this  in  the  boundary  condition  at  the  earth,  we  obtain 


_a 

3r 


v-  1/2  _  1/2 

<*7>  HV3  «i>  - k 


(tl> 


r  1/2  1/2 

(_1, 


-  r 


(6*32) 


r=a 


where 


k  =  -K  /K  . 
V1  V2 


(6*33) 


In  the  case  of  a  perfectly  conducting  earth  T 
(9-32)  reduces  to 


°°,  and  equation 


r  1/2  1/2 

(r  )  -  k  (— ) 
Sj/  1/3  U1; 


Hl/3 


=  0 


(6.34) 


r=a 


or 


H 


(2) 

1/3 


H 


(1) 

1/3 


<*!> 


r=a 


=  k  . 


(6.35) 


If  (6.22)  and  (6.23)  are  substituted  into  (4.14)  the  coefficient 
in  no  way  enter  the  expression  U(r,0). 
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SECTION  7 


MODIFICATION  OF  THE  RADIAL 
DIFFERENTIAL  EQUATION 

The  next  step  is  to  obtain  solutions  of  the  homogenous  differential 
equation  (3.5).  We  will  write  equ.  (3.5)  in  the  following  form: 


d"u  0  du  0  i \ 

n  ,  2  n  ,  r,  x  n(n+l)  t  n 

~T  *  r  57"  +  Ik  (r)  "  — - 1  un  '  0  ' 

dr  r 


l  /  •  1  f 


Letting  h  =  r-a  and  substituting  in  (7.1),  we  have 


d2u 


o-  +  [k  2  k2  (h)  -  — 1 5  ]  u  =  0 
dh2  °  (h+a)2  n 


2  2 
a  k 


V  7 . 2) 


2  2 

where  we  have  written  a  k  for  n(n+l) 

n 


Now 


2.2 

a  k 


(h+a) 


t  ,2  n  2hv 

■k  «  k  (1  -  — )  , 

2  n  a 


(7 .  Vi 


so  that  equation  (7.2)  becomes 


d2u 


dh 


n  r,  2  ,  2,,  .  .  2  ,  2  h.  n 

_ - 1-  [k  k  (h)  -  k  +  2k  — ]  u  =  0 

2  o  n  nan 


(7.4) 
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Since  usually, 


k(h)  «  1 


and 


1;  =  k 
n  o 


the  part  in  brackets  in  equation  (7.4)  becomes 


k  2  k2(h)  +  2k  2  -  =  k2N2(r) 

o  n  a 

2  2 

to  within  terms  of  the  order  of  h  /a  ,  and 

NCr)-^  • 

a 

Letting 

N2  (h)  =  1  +  y(h) 
y (h)  =  2  x  10“6  [M(h)  -  Ml 

and 


(1  -  xn) 


7 


we  obtain  for  equation  (7.1) 


d2u 


dh 


71  +  K  ty(h)  +  Xn]  un  =  0  * 


(7.5) 


(7.6) 


(7.7) 


(7.8) 


(7.9) 

(7.10) 


(7.11) 


(7.12) 
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The  effect  of  this  transformation  has  been  to  transform  the 
spherical  earth  problem  into  the  so  called  "flat  earth"  problem.  The 
resulting  differential  equation  (7*12)  turns  out  to  be  more  easily 
solvable.  In  the  case  when  the  index  of  refraction  varies  linearly 
with  height,  the  solutions  of  equation  (7.12)  take  on  particularly 
simple  forms.  The  price  one  pays  is  that  the  solutions  will  have  an 
error  of  approximately  2%  for  the  frequencies  and  elevations  of  interest 
here.  This  we  will  tolerate. 

In  terms  of  the  function  in  the  general  solution  (5.10),  equ. 
(7.12)  becomes 


d  Vh)  2 

- +  k  [y (h)  +  X  1  Q  (h) 

dh 


(7.13) 


The  general  solution  of  equ.  (7.13)  is 

1/2 


Vh>  ' 


[Av  (^h))  +  Bv  (ah))]  (7.14) 


where 


ah) 


(h  )  dh 


(7.15) 


and 


♦(h)  =  kQ  [y (h)  +  Av]  .  (7.16) 

Next,  we  will  write  the  Hankel  functions  of  order  one-third,  H^^(£) 
in  terms  of  the  Furry  modified  Hankel  functions  h(z) .  Some  of  the 
properties  of  these  functions  are  given  in  Appendix  E.  Thus,  the 
transformation  is 
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hi(z)  =  e1/3  (h)  (c(h)) 


(7.17) 


where 


2/3 


z  =  (-j  ?(h)) 


3/2 


C(h)  =  3-2 


and  equ.  (7.14)  becomes 

.1/6 


Vh)  -  T72Tr  [AV  hi  +  Bv  h2  (*» 

<P  (n) 


(7-18) 

(7-19) 

(7-20) 


This  form  of  the  solution  to  equation  (7.1)  will  be  used  in  all  sub¬ 
sequent  calculations.  The  remaining  equations  of  Section  5  and  Sec¬ 
tion  6  are  transformed  in  a  straight  forward  if  somewhat  tedious  way. 
In  particular  equation  (6.35)  now  takes  the  form 


h2(z(h)) 

h^(z(h)) 


h=0 


k  . 


(7.21) 


38 


SECTION  8 


MODEL  FOR  THE  TROPOSPHERE 


As  we  have  seen  previously,  a  smoothly  varying  index  of  refraction 
with  an  elevated  layer  will,  in  general,  result  in  the  occurrence  of 
three  transition  points.  Instead  of  an  analytic  function  for  the  index 
of  refraction,  we  will  choose  three  linear  line  segments  which  inter¬ 
sect  at  the  heights  h^  and  as  in  Figure  8.1.  The  heights  h^  and 

h^  correspond  to  those  of  n ^  and  in  Section  6.  These  are  the 
height  at  which  we  '’joined"  the  solutions  in  the  three  regions. 

is  the  height  of  the  layer  and  (h^  -  h^)  is  its  thickness. 
Mathematically,  the  model  can  be  described  as, 


y(h)  - 


+  aih 


0  1  h  1  hi 


y  L  +  a2(h  -  h1)  £  h  £ 


y 2  +  ^(h  “  h2)  h  1  h2  • 


(8.1) 


We  have  chosen  the  same  slope  (standard)  for  y(h)  in  region  I  and  III. 
The  slope  a2  in  region  II  is  simply  given  by 

y(h2)  -  y(h1) 


The  model  y(h)  is  illustrated  in  Figure  8.2. 

For  this  choice  of  model  the  variables  £  in  equ.  (6.4)  to  (6.6) 
and  (6.14)  to  (6.17)  take  on  the  form 
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Figure  8.1  LAYER  MODEL 


Figure  8.2  THE  FUNCTION  Y(h). 


r  h  1/2 

4n  (h)  =  k  /  [y  +  a.h  +  X  ]  dh 

1  of  Jo  1  n 

h01 

(8.3) 

r h  1/2 

=  +  a2  ^  +  A  ]  dh 

(8.4) 

rh  1/2 

C-(h)  =  k  /  [y0  +  a  (h  -  h0)  +  A  ]  dh  , 

i  o  /  z  I  2  n 

h03 

(8.5) 

and 


7  k  3/2 

r’ll(hl)  =  3  (a^  [yl  +  Xn^ 

(8.6) 

?  K  3 /2 

WV  "  I  ^  tyl  +  V 

(8.7) 

,  k  3/2 

*22  <V  =  T  ^2  +  V 

(8.8) 

C32<h2>  '  I  ‘y2  +  Xn]3/2' 

(8.9) 

In  terms  of  the  variable  z  defined  in  equ.  (7.18),  we  have  for  equ. 

(8.6) 

to  (8.9) 

k  2/3 

zn  ■  ^  <yi  +  V 

k  2/3 

Z21  ’  (yl  +  V 

(8.10) 

(8.11) 
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(8.12) 


k  2/3 

z22  =  <“>  <y2  +  Xn> 

k  2/3 

Z32  "  <y2  +  V 


(8.13) 


Since,  for  the  model  of  Figure  8.2,  is  negative,  we  will  choose  the 
principal  root,  that  is 


-i|i  k  2/3 


Z21  "  6 


3  i  o 


^  (^1  +  V 


(8.14) 


and 


-iii  k 


z22  ‘  e 


3  i  o 


a2'  <y2  +  V  • 


(8.15) 


In  terms  of  this  model,  the  solution  Q^(h),  equ.  (7.20),  takes  the 
following  form  in  region  I • 


1/6 


Qu00  -  t-i-jj 

3alko 


2/3 

\hi  (yc  +  aih  +  v) 


k  2/3 

+  Bvh2  (O  (yG  +  aih  +  v) 


(8.16) 


for  0  <_  h  h  ; 


in  region  II: 
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Vh>  =  w 

2  o 


k  2/3 

Cvhl  (yl  +  a2(h  -  V  +  v)  + 


(8.17) 


°vh2  (yl  +  a2(h  "  V  +  K) 


for  _<  h  <_  ; 


and  in  region  III  • 


1/6 


Vh>  ■  0 


3a1  k 
1  o 


F  h0 
v  2 


'k  2/3 

^  (y2  +  al(h  -  V  + 


(8.18) 


for  h  21  • 


Expressions  (8.16)  to  (8.18)  are  substituted  in  the  transcendental 
equation  (7.21)  and  solved  for  the  eigenvalues  A  . 
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SECTION  9 


EXPRESSION  FOR  THE  ELECTRIC  FIELD 


In  terms  of  the  function  and  variables  defined  in  Section  7  and  8, 
the  expression  for  the  potential  U(r,0)  in  equation  (4.13)  takes  the 
form 


U(h,0) 


2v  +  1 
sin  v0 


v 


Q  (h  )  Q  (h)  P  (cos(tt  -  0)) 
V  O  V  V 


P_d 

dh 


Qv(h) 


h=0 


_8 

at 


Vh> 


h=0 

t-v 


(9.1) 


where  v  is  related  to  X 

v 


by 


v 


koa(1  "  I  V 


(9.2) 


and 


_a 

at 


-2  a 


k  a  ax 
o  v 


(9.3) 


and  h  refers  to  the  height  of  the  transmitter, 
o 

The  coefficients  A  and  B  can  be  related  to  k  using  equ.  (6.22) 

v  v  & 

and  (6.23)  and  (6.32).  Thus  the  expression  for  the  height  function 
Q^(h),  equ.  (8.16)  becomes  in  region  I 


Q  (h)  -  F  K  ( - 5-) 

VVo  2 

z  3a,  k 
1  o 


1/ 6r  /k  2/3 

ho  ( (~)  (yQ  +  axh  +  xv)|  - 


(9. A) 


/k  2/3  \ 

k  hl  ^  (*o  +  +  V, 
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for 


0  1  h  1  hx  • 


Similarly,  the  coefficients  C  and  D  can  be  written  in  terms  of 

v  v 

F  K ^ •  However,  in  the  present  investigation  we  will  limit  our¬ 
selves  to  determining  the  fields  above  and  below  the  layer.  For 
h  _>  h2 ,  we  have 


Qv(h)  -  Fv  (■ 


3a,  k 
1  o 


1/6  /k  2/3  \ 

'  h2  (^2  +  al(h  "  h2)  +  V/ 


(9.5) 


for 


h  >  h2 


Substituting  equation  (9,4)  into  (9,1)  and  carrying  out  the  necessary 
differentiations,  we  obtain  the  following  expression  for  the  potential 
below  the  layer, 


-M 


,  2  4/3,,  x  1/3 
4a  a  (k  p) 
o 


sin  vtt 


P  (cos  (tt  -  0) 


(h1(0)h2(h0)  -  h2(0)h1(h0))(h1(0)h2(h)  -  h2(0)h1(h))/ 


1 


a  (k  ap  ) 
o 


2/3 


(9.6) 


for 


0  1  h  1  hx  • 
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The  notation  1^(0),  h^Ch^)  and  h^(h)  refers  to  the  value  of  the  Hankel 
function  at  the  values  of  h  =  0,  h  =  and  h  =  h,  respectively,  in  equ. 
(9.4).  In  this  derivation  use  was  made  of  the  Wronskian 

hl(?l)  h2f  (Cx)  -  (?x)  h2  C^)  =  -cti  (9.7) 

where  a  =  1.457495.  The  coefficient  p  is  related  to  a^  by 


a 


1 


1_ 

pa 


(9.8) 


p  is  equal  to  (2/3)  for  the  standard  (4/3)  earth  atmosphere.  In  equ.  (9.6) 
kT  is  the  derivative  of  k  with  respect  to  the  eigenvalue  A^. 

Since  v  is  of  the  order  of  k^a,  we  can  make  use  of  the  asymptotic 
expressions  for  the  Legendre  functions  of  large  order. 

In  that  case. 


(v  +  y)  (cos  (tt  -  0)) 
sin  vtt 


3  ( 


2k  a 
o 

tt  sin  6 


1/2 

)  exp 


.  3ir  .  ,  ,  1. 

1  —  -  1  (V  +  y) 


(9.9) 


In  future  discussions  it  will  be  of  interest  to  compare  the  potential 
U  to  its  free  space  value.  The  value  of  U  is 


w  -ik(b)R 

u  =  —  * - 

o  4tt  R 


(9.10) 


If  we  approximate  the  distance  R  by  a6  and  make  use  of  equ.  (9.9)  in  equ. 
(9.6),  we  obtain 
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i  /  /-j  i  (  3tt  .  1  « 

U_  _  .  9  N1/2  1  4~  1  2  6 

U  (sin  0^  e  6  2  1/3 

o  a  p 


2 it  (k  a)1/30 

O 


1/2 


E 


i  k  a  0  X 
2  o  v 


(h^ (0)h2 (hQ)  -  h2(0)h1(h0))(h1(0)h2(h)  -  h2(0)h1(h)) 


1  -  ~  ~  -.2/3  k’  "l2  (0) 

a(k  ap) 
o 


(9.11) 


for  0  <_  h  _<  h^.  The  eigenvalues  X  are  obtained  from 


-  Cj 

,  te.”0 


(9.12) 


where  k  is  given  by  equation  (6.32). 

The  electric  field  E,  is  obtained  from  U  by  using  the  expression 

<P 


.  3U 
%  =  -iUy  30 


(9.13) 


in  equ.  (2.2).  Since  0  appears  only  in  the  asymptotic  expression  (9.9), 
we  have 
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_3 

36 


1/2 

ir  sin  6  I  exP 


2k  a 
o 


if1-!  (v  +  j)  6 


(9.14) 


2k  a 
o 


V7T  sin  0 


\ 1/2 

.) 

-ik  a 
o 

exp 

3  ir  1 

i  f-  -  i  (V  +  f)  e 

- 

- 

Using  equ.  (9.14),  we  obtain  for  the  electric  field  E .  , 

<P 


Mko2/3uy 

E<()  =  ,_2_1/ 3_  1/ 3 


4a  a  p' 


2  1/2  3W  i 

^2ko  \  1  r  ~ 1 1 0 

,  7T  sin  0  /  e  e 


2 


i  k  a  0  X 
*  2  o  v 


(9.15) 


(hx (0) h2 (hQ)  -  h2(0)h1(hQ))  (h1(0)h2(h)  -  h2(0)h1(h)) 


1  - -  273  k’  hl2(0) 


a(k  ap) 
o 


for  0  <_  h  _<  . 

It  should  be  noted  that  the  ratio  U/Uq  does  not  differ  appreciably  from 
the  ratio  of  electric  fields  E./E.  .  As  long  as  the  distance  R  at 

which  the  field  is  measured  is  very  much  greater  than  a  wavelength  it 
is  sufficient  to  use  the  ratio  U/Uq. 
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It  will  be  of  interest  in  later  calculations  to  compare  the 
value  of  the  potential  U  in  a  layered  atmosphere  with  those  values 
when  no  layer  is  present.  For  an  homogenous  troposphere  we  then  have 


U 

U 


(sin  0) 


2tt  (k  a)^^0 
_ o _ 

1/3 


1/2 


3tt 

4 


(9.16) 


z 


i  tt  k  a  0 
2  o 


A  v  h2^ho^h2^h^ 


[h2'(°)J 


where  the  A are  roots  of  the  equation 


2/3 


h2  V  =  0 


and 


/k 2/3  \ 

h2<V-  h2  %  +  alh0  +  V) 


(9.17) 


(9.18) 


As  a  check  on  these  results,  it  is  of  interest  to  compare  the 
term  in  the  denominator  of  equation  (9.11),  namely 


1  -  — ‘ -271  k\2(V 

a  (k  ap  ) 
o 


(9.19) 


with  the  result  of  Furry  [16]  for  completely  trapped  modes.  Here  we 
have  written  h^(A^)  instead  of  h^(0)  to  emphasize  the  fact  that  A^  is 
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complex.  This  factor  is  proportional  to  the  normalization  constant 
for  the  eigenfunctions,  i.e.,  the  solutions  of  equation  (7.13).  For 
this  factor  Furry  gets,  in  terms  of  our  notation 


(9.20) 


where  the  eigenvalues 


X  are  written  as 

v 


X 

v 


a  +  ig 
v  v 


(9.21) 


Carrying  out  the  integration,  we  get 


-1/3 


-2 


LP  a 


KM 


1  - 


1/3  (iia 

°  V  ala2  / 

,  2  2,1/3 
a(p  a  ) 


(yi  +  %)1/2/hl(%)2 


(9.22) 


We  will  concentrate  on  the  factor  in  large  parenthesis.  For  completely 
trapped  modes  we  obtain 


k  —  e 


iir/3 


i 


3/2 

K  +  xv] 


(9.23) 
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so  that  expression  (9.19)  becomes 


1  - 


2h  1/3  (VJi k  +  X  1 

a  \  axa2  /  1  V 


1/2 


+  if 


a(ap) 


2/3 


(9.24) 


-  1  ^  fc— "A  (y.  +  x  )1/2  h  2(X  ) 

3  \  a1a2  J  1  »  1  “ 


But  at  the  eigenvalue  for  completely  trapped  modes  (see  Section  10.2) 


(\zr)  (yi +  V3/2  -  (m  - 1/2)*  m  =  1’2»3’ 


(9.25) 


Then  expression  (9.24)  becomes 

1/3  /a2  ‘  al\ 

2k  — - -  ( y.  +  X  ) 

°  \  ala2/  1  V 


1/2 


1  + 


a(ap) 


2/3 


V"  w* 


(9.26) 


For  the  completely  trapped  modes  A^  =  ,  to  a  high  approximation. 

Furthermore,  Furry  introduces  the  correction  factor  6  so  that 


iir/6_  i6 

e  hj(a^)  =  pure  imaginary 


(9.27) 


Using  this  fact  in  equation  (9.26)  we  obtain, 
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1/3 


2k 


1  - 


'a2  ~  V 

V  ala2/ 


(yl  + 


1/2 


a(ap) 


2/3 


!hn (a  ) 

1  1  v 


(9.28) 


which  agrees  with  equation  (9.22)  for  the  term  in  brackets. 
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SECTION  10 
NUMERICAL  RESULTS 

10 . 1  Description  of  the  Computer  Calculations 

The  computations  were  done  in  three  major  steps.  First,  the 
Hankel  functions  were  calculated;  second,  the  roots  of  the  transcen¬ 
dental  equation  (7.21)  were  found  and  finally  the  Hertz  potential 
function,  essentially  a  sum  of  the  roots,  was  obtained.  In  the  ex¬ 
pressions  for  the  roots  of  the  transcendental  equation  and  the  Hertz 
potential  or  electric  field,  Hankel  functions  of  the  first  and  second 
kind  of  order  one-third  and  complex  argument  z,  i.e.,  H^^^(z)>  and 
their  derivatives,  appear.  The  variable  z  takes  on  the  whole  range  of 
complex  values,  from  very  small  to  very  large.  To  evaluate  these  func¬ 
tions,  we  made  use  of  the  Furry  fl6]  infinite  series  and  asymptotic 
series  expressions  for  the  modified  Hankel  functions  of  order  one- 
third,  tsee  Appendix  D  and  Section  10.2].  As  our  condition 

on  the  accuracy  of  the  calculations,  we  required  agreement  with  the 
tabulated  results.  This  meant  that  the  modified  Hankel  functions  had 
to  be  calculated  at  least  in  double-precision  (fourteen  place  accuracy) 
on  the  IBM  370-155*.  (This  was  also  the  limitation  of  the  present 

software  available.)  It  turned  out  that  the  agreement  of  the  calcul¬ 
ated  results  with  the  tables  was  good  to  seven  (out  of  eight)  signifi¬ 
cant  figures  after  the  decimal  point.  This  was  satisfactory.  The 
series  expression  for  the  Hankel  functions  were  used  for  |z|  <6  and 
the  asymptotic  expressions  for  |z|  ^  6.  The  series  in  the  asymptotic 
expressions  was  carried  to  powers  in  z  of  z  .  The  agreement  between 

the  two  series  at  z  =  6  is  one-to-one.  In  the  case  of  the  asymptotic 
series  one  must  be  very  careful  so  as  to  use  the  proper  expressions  in 


As  an  historical  note,  it  is  interesting  to  reflect  that  the  tabulated 
results  were  obtained  by  the  Harvard  Computation  Laboratory  (now  Aiken 
Computational  Center)  in  1944. 
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the  various  domains  of  the  complex  plane.  In  particular,  if  the 
magnitude  of  z  is  large,  and  z  is  very  near  the  negative  real  axis, 
it  is  more  convenient  to  use  the  Miller  asymptotic  formulas  valid  for 
<  arg  <  y1-  [16]  . 


The  next  major  task  was  the  calculation  of  the  roots  of  the 
transcendental  equation  (7.21).  This  turned  out  to  be  the  most  diffi¬ 
cult  and  time  consuming  operation  of  the  whole  program.  The  problem 
was  complicated  not  only  by  the  fact  that  the  roots  were  in  the  complex 
plane  but  also  by  the  fact  that  roots  of  large  real  part  and  extremely 
small  imaginary  part  (trapped  modes)  had  to  be  obtained.  The  expression 
for  the  transcendental  equation  (7.21)  can  be  written  either  as 


F  =  h2(Av)  -  k(hrh2,Xv)  hx  (Av)  =  0 


(10.1) 


or  as 


W 

Fsorrk(Vh2>V  =  0 


(10.2) 


To  obtain  the  roots  corresponding  to  the  trapped  modes  it  is  essential 
to  use  the  second  form  for  F,  since  both  h2  and  h^  are  very  large 
simultaneously.  The  secant  and  the  Newton-Raphson  iteration  methods 
[17]  were  used  to  solve  F  =  0  for  A  .  The  experience  of  previous 
investigators  attempting  to  solve  similar  problems  had  been  that  the 
two  methods  give,  in  general,  good  convergence  in  the  complex  plane. 
Each  method,  however,  was  cited  to  give  no  convergence  for  particular 
geometries  of  F  in  the  complex  plane.  The  recommendation  in  those 
instances  has  been  to  go  simply  to  another  method.  It  would  be  pre¬ 
sumptuous  to  make  a  general  conclusion  on  which  of  the  two  methods 
gives  more  rapid  convergence,  for  this  depends  on  the  particular  pro¬ 
blems  involved;  speaking  only  for  our  problem  it  was  found  that  the 
Newton-Rhapson  method  gives  somewhat  better  rate  of  convergence.  The 


54 


roots  found,  however,  by  each  method  are  identical  for  the  same  error 
bound  on  A^ ,  i.e.,  |  A^n^  -  A^n  "^  |  <  e.  The  disadvantage  of  both 
methods  is  the  requirement  that  the  initial  guess  of  a  root  must  be 
fairly  close  to  the  actual  location  of  that  root  to  obtain  convergence 
at  all.  The  advantage  of  the  secant  method  is  that  it  involves  only 
F  whereas  the  Newton-Rhapson  method  involves  the  derivative  of  F,  which 
in  our  problem  is  not  a  trivial  calculation.  Since  we  were  guided  by 
minimizing  the  computer  running  time,  all  roots  presented  in  this  re¬ 
port  were  found  by  the  Newton-Rhapson  method. 

The  calculation  of  the  Hertz  potential  presented  no  special  com¬ 
putational  problems,  and  hence  needs  no  further  discussion.  A  more 
detailed  description  of  the  computational  work  described  above,  in¬ 
cluding  computer  flow  charts,  will  appear  in  a  companion  report, 

10. 2  Calculation  of  the  Characteristic  Values  for  the  Propagation 

Modes 

The  complex  roots  or  characteristic  values  of  the  transcendental 
equation  (7.21)  can  be  thought  of  as  propagation  modes  similar  to  those 
encountered  in  waveguide  theory.  However,  the  boundary  value  problem 
of  the  elevated  duct  is  much  more  complicated  than  that  since  the  ele¬ 
vated  layer  can  not  be  simply  modeled  as  a  perfectly  reflecting  stratum. 
The  modes  that  actually  exist  can  be  classified  as  "trapped”  modes, 
"transitional"  modes  and  "leaky"  modes.  This  becomes  clearer  if  we 
think  of  the  geometric  interpretation  in  Figure  10.1.  The  first  trapped 
mode  is  m  =  1  and  the  last  is  m  =  5 .  There  are  always  a  finite  number 
of  these.  In  terms  of  the  roots  A^  of  equ.  (7.21), 


vv 

iOTT  -  k(V  ’  0 


(7.21) 


This  means  that  the  imaginary  part  of  X^  is  extremely  small  and  to  a 
high  approximation  X^  =  -y(h^).  Mode  number  m  =  5,  for  example,  is 
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MA-37,817 
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Figure  10.1  PROPAGATION  MODES  (EXAGGERATED  DIMENSIONS) 


confined  essentially  to  the  heights  between  h^  and  h^  and  is  char¬ 
acterized  by  a  very  small  variation  of  attenuation  with  range.  When  the  imaginary 
part  of  A^  ,  Im(A^),  is  small  and  the  real  part  Re(X^)  ~  -yO^),  we  are 
said  to  have  transitional  modes,  e.g.,  m  =  6  in  Figure  10.1.  Physically 
these  modes  correspond  to  waves  that  just  leak  out  of  the  duct.  The 
characteristic  values  A^ ,  that  have  real  parts  comparable  in  magnitude 
to  their  imaginary  parts  are  said  to  be  leaky.  These  modes  do  not 
allow  such  a  simple  geometric  interpretation  as  the  previous  two. 


Now  let  us  turn  to  the  solution  of  equation  (7.21)  for  A^ .  We 
started  by  finding  the  roots  corresponding  to  the  trapped  modes.  But, 
as  we  pointed  out  in  Section  10.1,  for  rapid  convergence  to  a  root  by 
either  the  secant  or  the  Newton-Rhapson  iteration  methods,  it  is 
essential  that  the  initial  guess  of  a  root  must  be  fairly  close  to  the 
root  itself.  To  get  a  handle  on  this  problem  we  approximated  the 
Hankel  functions  appearing  in  (7.21)  by  their  asymptotic  representa¬ 
tions,  keeping  only  the  first  terms.  In  this  case  k(A^)  is  given 
approximately  by 


k  ^ 


4k  a?  -  a- 

-1  -f  <-^>‘*1  +  V 


3/2 


where  the  notation  has  been  explained  previously.  And  equation 
becomes  essentially 


(10.3) 

(7.21) 


1 


+  e 


-iX 


=  0 


(10.4) 


where  X  is  given  by 
4k 


a?  “  ai 

x  -  -T  (^zr->  [5,1 +  V 


3/2 


Equation  (10.4)  is  satisfied  for 


(10.5) 
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X  =  nnr 


(m  is  odd). 


(10.6) 


As  a  check,  we  compared  this  result  with  the  Furry-Gamow  phase- 
integral  method  [  4  ] .  In  this  method  the  real  part  of  X is  deter¬ 
mined  by 

h°2  1/2  . 

Iy(h)  +  Xv]  dh  =  (m  -  f)  tt  (10.7) 

Tn  —  1,2,3*  •  • . 

where  the  h^  and  h^  satisfy  y(h)  +  X^  =  0.  A  straight-forward 
integration  in  two  parts  from  h^  to  and  h^  to  h^  gives 


2k 
_ a 

3 


,a2  - 


ala2 


3/2 

[yl  +  V 


m 


(10.8) 


1,2,3. . . 


It  is  seen  that  the  results  are  identical.  Equation  (10.7)  is  known 
as  the  "phase-integral  condition".  The  cases  for  which  this  method 
applies  are  as  follows: 


(i)  The  curve  y(h)  must  have  at  least  one  minimum,  in  our 
case  it  has  two,  y(h^)  and  y(h^)>  a  mode  to  which  the 
method  can  be  applied  must  have 

-yG^)  <  Re(Av)  <  -y(h2>  (10.9) 

so  that  the  waves  are  trapped;  and  it  must  have 


Im(Xv)  «  -(y(h1)  +  Re(A^)) 
Im(Xv)  «  -(y(h2)  +  Re(A^))  ; 
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(ii)  The  value  of  Re(A^)  is  found  from  equ.  (10.8). 

After  the  Re(A^)  is  known,  Im(A^)  is  calculated 
from  an  explicit  real  expression  to  be  derived 
presently . 

For  the  determination  of  the  imaginary  part  of  A^  we  will  use 

the  modified  index  of  refraction  profile  y(h)  of  equ.  (8.1),  Figure 

8.2.  Let  us  write  the  characteristic  values  A  as 

v 


A  =  a  +  i(3 
v  v  v 


(10.11) 


We  will  now  extend  the  formulas  Furry  [3]  derived  for  a  single  minimum 
of  y(h)  to  two  minima.  The  imaginary  part  6^  is  then  given  by  (see 
Figure  10.2) 

0  -2  W. 

(10.12) 


o  .  -2  W 

e  =  C  2  k  1  e  v 
v  v  a 


where 


-2 


2f 


02 


1/2 


Cv  ■  k  2  |  (y(h)  +  a^)  dh 


W  =  k 
m  a 


01 


/ 


03  1/2 

I  y  (h)  -  aj  dh 


(10.13) 


(10.14) 


02 


and  k  *  2ir/A,  where  A  is  the  wavelength.  When  equs .  (10.18)  and 
a 

(10.14)  are  applied  to  the  profile  of  Figure  10.2  we  obtain 

-2  1/2  a2  "  a1 

C  =  4  (y.  +  a  )  (— - -) 

v  1  v'  a\a2 


(10.15) 
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and 


2k 


W  = 
m 


<-y2  -  «v) 


3/2 


a2  “  al 
(— - -) 

ala2 


The  are  determined  from  equ.  (10.6)  or  (10.8) 


we  shall  call  these  a 


8^  the  Furry-Gamow  roots, 


equ.  (10.16)  that  the  condition  (10.10)  requires 


(10.16) 


For  later  comparison 
It  is  seen  from 


and 


W 

m 


<< 


1 


«  1 


(10.17) 


This  implies  that  the  values  of  8  will  be  in  error  near  a  =  -y„  and 

v  v  J 1 

a  =  -y  .  It  is  surprising,  however,  how  good  these  approximate  values 

of  a  and  8  are  to  the  exact  values  determined  by  the  Newton-Rhapson 

v  v 

method.  As  a  particular  example  we  have  chosen  the  following  constraints 


X 

1  meter 

-6 

II 

211.66  x 

10  D 

-6 

y2  " 

131.66  x 

10 

-7 

al  = 

2.3518  x 

10  7 

■  7 

a2  = 

-8.0  x  10 

In  Table  10.1  are  shown  the  Furry-Gamow  roots  and  their  exact  values. 
From  these  results  it  is  seen  that  the  real  part  of  X is  in  error  by 
less  than  one  percent.  The  imaginary  part  of  X^  is  in  error  by  about 
10%;  except  for  m  =  2  where  the  error  is  about  30%.  The  10%  error  can 
be  explained  from  the  fact  that  condition  (10.17)  is  not  always  satis¬ 
fied.  There  is  no  ready  explanation  for  the  30%  error  in  mode  number  2. 
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TABLE  10.1  COMPARISON  OF  FURRY -GAMOW  ROOTS  AND  EXACT  ROOTS 


It  can  be  concluded,  however,  that  the  Furry-Gamow  values  serve  as  an 
excellent  guide  for  the  initial  guess  of  a  root  and  for  this  reason 
alone  are  indispensible . 


There  are  no  such  relatively  simple  formulas  that  give  the  roots 
corresponding  to  the  transitional  and  leaky  modes.  It  was  fortuitous 
that  the  location  of  these  roots  could  be  predicted  reasonably  accur¬ 
ately  by  linear  extrapolation  from  the  trapped  mode  root  values.  How¬ 
ever,  in  only  one  instance  could  the  location  of  three  successive  roots 
be  predicted  by  linear  extrapolation.  The  reason  for  this  failure  can 

be  seen  from  the  non-linear  character  of  6  as  a  function  of  mode 

v 

number  in  Figure  10.3.  It  is  perhaps  advantages  to  use  some  curve 
fitting  routine  to  predict  more  root  values  in  advance;  this  was  not 
tried  here.  Such  a  procedure  would  certainly  cut  down  on  the  computer 
running  time  if  it  proves  successful. 


There  was  another  way  in  which  we  went  about  locating  a  root 
for  a  particular  set  of  values  of  the  parameters  A,  y^ ,  y^,  a^  and  a^ • 
This  approach  was  motivated  by  Northover's  [18]  investigation  in 
which  he  considered  reflections  from  an  elevated  layer  by  treating 
such  a  layer  as  a  discontinuity  in  the  index  of  refraction.  He  pre¬ 
dicted  that  for  the  important  modes  in  ground  to  ground  communications, 
namely  those  waves  that  are  reflected  by  the  elevated  layer  back  to 
the  receiver  should  have  corresponding  characteristic  values  whose 
magnitude  is  less  than  unity.  In  terms  of  our  notation  this  condition 
means  that 


<  1  • 


(10.18) 


They  correspond  to  leaky  modes.  It  turned  out  that  these  roots  could 
be  found  rather  quickly  since  there  were  always  just  a  few  of  them 
(see  e.g.,  Table  10.2).  And  the  real  and  imaginary  parts  were  of  com¬ 
parable  magnitude.  These  roots  then  provided  another  starting  point 
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CHARACTERISTIC  VALUES  AS  A  FUNCTION  OF  MODE  NUMBER 


for  finding  additional  characteristic  values.  This  procedure  proved 
to  be  very  convenient. 

10. 3  Discussion  of  the  Characteristic  Values  and  Their  Physical 

Interpretation 

All  the  characteristic  values,  X  ,  we  have  found  are  listed  in 

v 

Tables  10.2,  10.3  and  10.4.  The  accuracy  of  these  values  obtained  on 
the  computer  was  actually  14  places;  they  have  been  rounded  off  to  five 
places  after  the  decimal  in  the  tables.  The  values  were  obtained  as 
a  function  of  frequency  (wavelength  X),  and  layer  height,  hr,  as  a 
function  of  lapse  rate  within  the  layer,  and  outside  the  layer, 
a^  and  as  a  function  of  the  value  of  y(h)  at  the  start  of  the  layer, 
y^  and  at  its  termination  y^ . 

The  real  part  of  the  characteristic  values  enters  as  a  phase 
factor  in  the  final  expression  for  the  field,  equ.  (9.15).  Each 
characteristic  value  or  mode  is  seen  to  have  a  different  phase  and 
hence,  the  addition  of  many  such  modes  results  in  appropriate  phase 
vector  addition.  The  imaginary  part  of  X  corresponds  to  the  attenua¬ 
tion  of  the  mode.  This  can  be  written  as 

1/2  6  k  d 
v  o 

attenuation  =  — - 7— -  (10.19) 

It  is  seen  from  the  tables  that  the  trapped  modes  have  a  very  small 
imaginary  part  so  that  the  attenuation  of  these  is  like  (l/^^d).  It 
is  noted  from  a  comparison  of  Tables  10.2  and  10.3  that  the  thicker 
layer  (100  meters)  has  the  smaller  imaginary  part  for  the  first  mode, 
hence  the  attenuation  of  this  mode  is  smaller.  It  will  be  seen  in 
Section  10.4  that  this  will  also  correspond  to  more  energy  being 
trapped  in  the  duct.  In  both  cases  the  first  five  modes  are  trapped. 

In  Table  10.4,  the  layer  is  even  thicker  (150  m)  but  we  have  lowered 
the  frequency  to  150  MHz  and  as  a  result  the  first  mode  has  an  imaginary 
part  several  orders  of  magnitude  greater  than  for  300  MHz.  In  addition, 
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Mode 

Number 

Characteristic 

Values 

m 

Real  Part  (x  10'5) 

Imaginary  Part 

i 

-19.25747 

.31597  x  10_13 

2 

-17.33864 

.56905  x  10 "10 

3 

-15.78268 

.77958  x  10-8 

4 

-14.43457 

.17336  x  10~6 

5 

-13.20194 

.84462  x  10-6 

6 

-12.00625 

.16483  x  10” 

7 

-10.85324 

.22246  x  10-5 

8 

-  9.75080 

.26214  x  10“5 

9 

-  8.69502 

.29129  x  10-5 

10 

-  7.68003 

.31404  x  10-5 

11 

-  6.70058 

.33267  x  10-5 

12 

-  5.75236 

.34849  x  10-5 

13 

-  4.83184 

.36233  x  10“5 

14 

-  3.93611 

.37476  x  10-5 

15 

-  3.06264 

.38639  x  10-5 

16 

-  2.21008 

.42040  x  10-5 

17 

-  1.36162 

.42117  x  10-5 

18 

-  .49410 

.46243  x  10” 

19 

.43138 

.52265  x  10-5 

20 

1.43287 

.59324  x  10 

21 

2.51086 

.66949  x  10“5 

22 

3.66053 

.75062  x  10“5 

23 

4.87641 

.83685  x  10-5 

24 

6.15302 

.92792  x  10~5 

25 

7.48485 

1.02197  x  10-5 

26 

8.86721 

1.11422  x  10-5 

27 

10.29907 

1.19673  x  10-5 

28 

11.78600 

1.26204  x  10-5 

29 

13.33784 

1.30910  x  10-5 

30 

14.96139 

1.34344  x  10-5 

TABLE  10.2  CHARACTERISTIC  VALUES 
(A  =  1m,  AM  =  40  in  30m,  =  900m) 
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Mode 

Number 

1 

Characteristic  Values 

I 

m 

Real  Part  (x  10  5) 

'i 

Imaginary  Part 

i 

-19.58371 

.12851  x  10"15 

2 

-17.74367 

.77288  x  10-12 

3 

-16.28801 

.32337  x  10-9 

4 

-15.06395 

.22680  x  10-7 

5 

-14.00183 

.32445  x  10"6 

6 

-12.98671 

.12363  x  10"5 

7 

-11.94416 

.22376  x  10-5 

8 

-10.90414 

.30204  x  10"5 

9 

-  9.89299 

.36548  x  10“5 

10 

-  8.91984 

.41982  x  10-5 

11 

-  7.98680 

.46395  x  10"5 

12 

-  7.08954 

.49208  x  10"5 

13 

-  6.21641 

.50117  x  10"5 

14 

-  5.35600 

.49691  x  10-5 

15 

-  4.50455 

.48825  x  10-5 

16 

-  3.66337 

.48082  x  10-5 

17 

-  2.83430 

.47709  x  10“ 5 

18 

-  2.01671 

.48026  x  10“5 

19 

-  1.20037 

.49915  x  10"5 

20 

-  .35700 

.54292  x  10"5 

21 

.54346 

.60849  x  10"5 

22 

1.51098 

. 68840  x  10-5 

23 

2.53954 

.77999  x  10-5 

24 

3.61628 

.88048  x  10"5 

25 

4.72710 

.97651  x  10-5 

26 

5.87686 

1.04212  x  10"5 

27 

7.09619 

1.07111  x  10"5 

28 

8.40164 

1.08675  x  10-5 

TABLE  10.3  CHARACTERISTIC  VALUES 
(A  =  lm,  AM  =  40  in  100m,  =  900m) 
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only  the  first  three  modes  are  trapped  resulting  in  lower  signal 
strength  within  the  duct. 

10.4  Comparison  of  Calculated  and  Experimental  Results 

The  calculations  of  field  intensity  relative  to  free-space  were 
performed  using  expressions  (9.11)  and  (9.12).  The  results  are  plotted 
in  Figures  10.4  to  10.8.  All  curves  are  for  a  frequency  of  300  MHz, 
and  a  layer  height  of  900  m.  Two  values  of  layer  thickness,  30  m  and 
100  m,  but  of  the  same  decrement  in  modified  index  of  refraction, 

AM  =  40,  have  been  used.  Such  layer  values  are  representative  of  those 
measured  in  the  trade  winds  ducts  near  San  Diego,  Calif.  [19]  and  the 
South  Atlantic  [201  .  The  field  intensity  is  plotted  as  a  function  of 
height  and  distance  for  various  transmitter  locations.  A  comparison 
of  the  plots  shows  that  the  thicker  layer  (100  meters)  has  also  the 
greater  field  intensity  values;  a  difference  of  about  10  dB  on  the 
average  exists  for  the  two  layer  thickness  values.  This  appears  logi¬ 
cal  since,  for  the  same  decrement  (AM  =  40)  the  100  meter  layer  results 
in  a  wider  duct  than  the  30  meter  layer  and,  hence,  is  capable  of 
trapping  more  energy.  When  the  transmitter  is  moved  out  of  the  duct, 
Figure  10.6,  the  field  intensity  decreases  by  approximately  5  dB  on 
the  average.  In  addition,  we  notice  in  this  case  interference  between 
the  waves  that  move  above  and  are  coupled  into  the  layer  and  those  that 
go  through  the  duct  strike  the  earth  and  then  escape  into  the  space 
above  the  layer.  In  Figure  10.6,  constructive  and  destructive  inter¬ 
ference  are  observed  at  approximately  450  km  and  500  km,  respectively. 

The  field  intensity  is  given  as  a  function  of  height  in  Figures  10.7 
and  10.8.  Here  again,  the  wider  duct  results  in  a  higher  field  strength, 
about  10  dB.  We  also  notice  in  Figure  10.7  that  we  obtain  interference 
type  (multipath)  lobing;  the  intensity  of  which  decreases  as  one  gets 
closer  to  the  bottom  of  the  layer.  Furthermore,  the  field  strength  in¬ 
creases  by  approximately  5  dB  from  near  the  bottom  of  the  duct  to  the 
top.  This  latter  behavior  does  not  appear  to  be  in  evidence  for  the 
wider  duct  in  Figure  10.8. 
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Figure  10.4  FIELD  INTENSITY  AS  A  FUNCTION  OF  TERMINAL 
SEPARATION  {  AM=  40  IN  30m  ) 


=  40  UNITS  IN  100  m 
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SEPARATION  (AM  =  40INI00m 
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Fig  u  re  10.6  FIELD  INTENSITY  AS  A  FUNCTION  OF  TERMINAL 
SEPARATION  (  A M  =  40  IN  1 00  m  ) 
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Figu  re  10.8.  FIELD  INTENSITY  AS  A  FUNCTION  OF  HEIGHT 

(  AM  =40  IN  100m) 
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We  will  now  compare  these  results  with  some  experimental  evidence. 

For  comparison,  we  use  the  test  results  of  Project  Neptune  [21].  This 
project  involved  air-to-air  and  air-to-ground  propagation  testing  in 
the  trade  wind  ducts  at  220  MHz.  Its  purpose  was  to  explore  the  communi¬ 
cation  potentialities  of  the  trade  winds  inversion  layers.  For  compari¬ 
son  we  will  choose  the  thinner  layer  (30  m) .  It  is  to  be  noted  that  in 
these  test  results,  the  measured  values  of  refractive  index  values  are 
open  to  question  [21  |.  The  agreement  between  our  results  and  the  tests 
is  good  on  a  relative  basis.  When  the  results  are  compared  on  an  absolute 
basis  our  results  are  approximately  10  -  15  dB  higher.  To  partially  ex¬ 
plain  this  difference,  we  note  that  our  results  are  for  a  somewhat  higher 
frequency,  300  MHz  versus  220  MHz.  In  general,  the  higher  the  frequency, 
the  greater  the  extent  of  trapping,  and  hence,  the  higher  signal  level. 
Another  difference  is  due  to  the  fact  that  in  our  model  the  duct  width 
and  AM  scenarios  constant  whereas  in  the  tests  the  width  and  decrement 
actually  changed  in  the  distance.  This  would  also  increase  the  leakage 
from  the  duct  and  result  in  a  lower  signal  level.  In  the  tests  it  was 
observed  that  signal  level  increased  from  5  to  20  dB  as  the  receiver 
moved  closer  to  the  bottom  of  the  layer.  The  field  being  highest  at 
the  base  of  the  inversion.  These  results  agree  essentially  with  ours. 

We  obtain  an  increase  in  signal  level  of  5  to  10  dB.  Furthermore,  it 
was  reported  that  the  signal  level  remained  "relatively  high"  above  the 
duct.  This  also  agrees  with  our  calculated  results.  It  being  that  the 
signal  is  only  approximately  5  dB  less  above  the  duct.  Finally,  the 
tests  established  good  signal  correlation  for  60  feet  ( ss  18  m)  verti¬ 
cally.  If  this  result  is  compared  with  Figure  10.7,  it  is  seen  that 
it  agrees  quite  well  with  the  calculated  numbers;  our  results  showing 
about  80  feet  at  the  3  dB  points. 
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SECTION  11 


SUMMARY  OF  CURRENT  RESULTS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 


We  will  here  summarize  the  results  we  have  obtained  and  explain 
their  significance.  Also,  we  will  make  several  recommendations  for 
future  work.  In  addition,  we  will  point  out  the  most  significant 
elements  to  be  incorporated  in  a  propagation  experiment. 

We  have  shown  that  the  presence  of  elevated  ducts  can  produce 
significantly  higher  signal  strengths;  on  the  order  of  20  dB  above 
free  space.  We  have  further  demonstrated  that  the  tri-linear  model 
for  the  index  of  refraction  is  plausible  and  gives  good  agreement  with 
experimental  results  on  a  relative  basis  but  tends  to  make  rather 
high  estimates  of  signal  strength  on  an  absolute  scale.  We  have 
shown  the  existence  of  interference  type  lobing  at  the  duct  height 
with  5  dB  to  10  dB  fades  at  300  MHz  and  good  signal  correlation  for 
vertical  distances  of  roughly  80  feet.  From  an  operational  point  of 
view,  it  is  best  to  fly  the  transmitter  and  receiver  at  the  base  of 
the  elevated  layer.  At  this  location  signal  strength  is  highest 
and  interference  type  lobing  is  least.  In  order  to  avoid  the  10  dB 
fades  it  is  essential  that  space  diversity  reception  be  employed. 

In  the  present  report  we  have  limited  our  analysis  to  the  detailed 
analysis  of  a  single  frequency,  300  MHz.  Future  work  should  be  ex¬ 
tended  to  a  comprehensive  comparison  of  several  frequencies.  Spatial 
correlation  distances  of  the  signals  should  be  obtained  for  such  fre¬ 
quencies  in  order  to  obtain  estimates  of  vertical  and  horizontal 
antenna  diversity  distances.  Furthermore,  a  spectral  analysis  should 
be  performed  to  estimate  the  signalling  bandwidth  that  can  be  supported 
in  the  presence  of  elevated  layers.  And  finally,  the  calculated  re¬ 
sults  that  have  been  obtained  for  heights  within  the  duct  and  slightly 
above  it  should  be  extended  to  include  the  whole  range  from  the  ground 
to  great  heights. 
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From  the  above  discussion  it  becomes  quite  clear  what  the  most 
essential  elements  of  any  air-to-air  propagation  experiment  should  be. 

The  most  detrimental  effect  of  this  channel  is  interference  type 
lobing.  Hence,  the  experiment  must,  first,  establish  the  signalling 
bandwidth  it  can  support.  This  could  be  performed  by  probing  the  channel 
with  a  multi-tone  signal  consisting  of  a  group  of  equally  spaced  fre¬ 
quencies.  Secondly,  the  experiment  must  test  some  antenna  diversity 
system  in  order  to  reduce  the  amplitude  modulation  caused  by  fading. 

A  general  experiment  incorporating  these  elements  is  described  in 
Reference  1. 
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APPENDIX  A 


JUSTIFICATION  FOR  THE  HERTZ  POTENTIAL 


We  shall  show  here  that  the  one  dimensional  Hertz  vector  n  = 

(0,0, U)  is  a  satisfactory  approximation  to  the  solution  of  Maxwell’s 
equations . 

Let  us  consider  II  which  satisfies  the  equation 
—  —  2  — 

v*v#n-vv*n-yeo)  n  =  o  (a-i) 

and  construct  the  field  vectors  corresponding  to  the  magnetic  type 
(magnetic  dipole) ,  then 


D  =  -  iy  cujV  •  n  (A-2) 

H  =  V  *  V  *  TT  (A-3) 

The  vectors  D  and  H  satisfy  the  equations 


V  x  E  = 


-  iioy  H  +  —  x  E 


V  x  H  =  iioe  E 


V  •  D  =  0 


V  *  B  =  Vy  *  H 


(A-4) 
(A-5) 
(A- 6) 
(A- 7) 


It  is  seen  that  equations  (A-4)  to  (A-7)  do  not  satisfy  Maxwell's 
equations  exactly.  For  the  atmosphere,  however,  the  gradient  of  y 
is  extremely  small  so  that  terms  involving  it  can  be  neglected.  The 
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ratio  of  the  magnitude  of  the  second  term  to  the  first  term  on  the 
right  hand  side  of  equ.  (A-A)  can  be  written  as 


X_ 

2tt 


1  vy I  .  yVo~  |e| 


(A- 8) 


Except  quite  close  to  the  source,  the  field  will  be  almost  a  plane 
wave  for  which  |e|/|h|  is  approximately  ^  y/e .  Furthermore,  |Vy|/y 
is  of  the  order  of  10“^,  and  therefore  the  ratio  of  the  second  term 
to  the  first  term  is  of  the  order  X  •  10“-*  ,  which  is  always  very 
small  compared  to  unity.  From  equs.  (A-5)  and  (A-7)  it  follows  that 

v  •  B  J_VyJ_  v/Vo  HL  /a_q \ 

2.  p  e  |i|  W  9) 

and  this  again  is  of  the  order  of  X  *  10  It  follows  that  equs.  (A-4) 

to  (A-7)  differ  from  MAxwell's  equations  by  a  negligible  amount,  and 
we  may  regard  equations  (A-2)  and  (A-3)  as  a  valid  way  to  obtain  the 
electromagnetic  field. 

Because  the  Hertz  vector  has  only  its  radial  component  non-zero, 
it  follows  that  [19] 


|v  •  b| 


V(ye)  •  V  •  II  =  0  and  V(ye)  *  II  —  0  •  (A-10) 
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APPENDIX  B 


CONTOUR  INTEGRAL  REPRESENTATION  FOR  U(r,0) 


Following  a  method  due  to  Watson,  we  shall  show  that  the  series 
expression  for  U(r,0),  (3.9)  can  be  expressed  as  a  contour  integral. 
Consider  the  series 


S 


(2n  +  1)  f (n)  P  (cos  0) 
n 


(B.l) 


where  f(n)  is  an  analytic  function  of  n.  We  will  show  that  S  can  be 
represented  by  the  following  contour  integral: 


■h 


—  f  (s  -  -r)  P  .  (cos  (tt  -  6)) 

COS  TTS  Z  s  -  1/2 


(B .  2) 


where  C 


C^  +  C^  is  a  contour  that  starts  at  00  -  i6  in  the  complex 


s-plane,  goes  below  the  real  axis  at  s  =  0  and  then  above  the  real  axis 
to  00  +  i6.  Since  f  (s  -  y)  and  Pg  _  (cos  (^  ~  8))  are  analytic  functions 

of  s,  the  only  singularities  of  the  integral  are  poles  at  those  values  of 
s  inside  C  for  which  cos  tts  =  0,  that  is  s  =  1/2,  3/2,  5/2,  ... 

Equation  (B.2)  can  be  written  as  the  sura  over  C^  and  C^ .  Hence, 


(B .  3) 


where  H(s)  is  the  integrand  in  (B.2).  In  the  second  integral  of  (B.3) 
replace  s  by  -s.  This  transformation  reflects  the  contour  C^  about  the 
origin  into  the  contour  C^,  as  well  as  alters  the  integrand. 
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Figure  B-l.  COMPLEX  S-PLANE 
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Hence,  equation  (B.3)  becomes 


H(s)ds 


to  (B.4)  results  in 


H(s)ds 


[H(s)  +  H(-s)]  ds  • 


Using  the  fact  that 


(B.4) 


(B.5) 


ps  _  i/2^cos  ^  =  P_s  _  ^  “  ®)) 


(B.  6) 


and  substituting  for  H(s)  we  obtain 


/  , 

JC1  +  S 


Ms‘(irs)  f  (S  ‘  I}  Ps  -  1/2 (coS  (1T  "  6))  ds 


+[  i  cos  to  [f  (_S  "  ¥  -  Hs  -  b]  ds 


(B .  7) 


The  contour  may  be  closed  by  a  semicircle  of  infinite  radius  in 

the  lower  half  plane  as  shown  in  Figure  (B.2).  Thus  (B.7)  becomes 
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(B.  8) 


/ 

Jc.  4. 


c+c.+c. 

13  4 


i  COS  7TS 


f(s  -  -|o  Pg  _  (cos  (tt  -  e)) 


ds  +  i. 


where  is  the  second  integral  in  (B.7). 

Now  substituting  the  expression  for  f  obtained  in  the  previous 
section, 


f(s  -  ±)  = 


-M 


2  4ttW 


s  -  1/2 


%  -  1/2 (b>  Vs  -  1/2 <r)  r±b 


Vs  -  l/2(b)  Qs  -  l/2(r)  rib 


(B.9) 


where  Qg  _  .^(r)  and  Vs  -  1/2^  rePresent  analytic  functions  along 
the  positive  real  axis.  Substituting  (5.9)  into  the  expression  for 
U(r,0)  we  obtain 


u(r>e>=i?r 


/ 


sds  s  -  1/2 


COS  7TS  W 


Cl+C3+C4 


s  -  1/2 


Qs  -  l/2<b>  Vs  -  l/2(r)  rib 


Vs  -  l/2(b)  Qs  -  l/2<r)  r^b' 


(B. 10) 


And  if  we  replace  s  by  v  +  —  and  apply  the  Caudy  residue  theorem,  then 
(B.10)  becomes 


w  (v  +  ^)  P  (cos  (-TT  -  6)) 

^5 - — 


<IT> 


s=v 


Vb>  Vr) 


VvCh)  Qv(r) 


+  I, 


(B.ll) 
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where  v  represents  the  zeros  of  Wg  in  the  lower  half  of  the  complex 
s-plane . 

It  can  be  shown  that  if  the  earth  is  assumed  to  be  a  perfect  con¬ 
ductor,  then  the  integral  vanishes  identically.  If  the  earth  is  not 
a  perfect  conductor,  1^  will  not  be  zero.  However,  it  can  be  shown 
that  1^  is  small  compared  to  the  first  term  of  the  series  in  (B.ll) 
as  long  as  the  dielectric  constant  of  the  earth  is  large  compared  to 
the  dielectric  constant  of  the  atmosphere.  So  that  finally. 


U(r,0) 


#  E 


sin  Vi: 


(cos  (tt  -  0)) 
dW 

‘ar> 

s=v 


Qv(b)  VvCr> 


Q  (r)  V  (b) 

V  v 


r  <  b 


r  >  b 


(B.12) 
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APPENDIX  C 


DERIVATION  OF  EQUATIONS  (4.12)  and  (4.13) 


From  equation  (4.4),  consider  the  following  two  equations: 


and 


(rVt)  + 


t(t  +  1) 

2 

r 


(rVt)  =  0 


<rV  + 


v(v  +  1) 

2 

r 


(rQv)  =  0 


(C.1) 


(C.2) 


Here  t  is  any  value  and  is  not  an  eigenvalue.  Multiply  (C.l)  by  (rQ^) 
and  (C.2)  by  (rV  ) ,  subtract  and  integrate  with  respect  to  r  from  r  =  a 
to  r  =  °°.  We  obtain 


(v  -  t)  (t+v+1) 


>  /  Wr  ’ 


rV  (rQ  ) 

t  dr  v  xv 


(rQ  )  -r-  (rV  ) 
xv  dr  t 


(C.  3) 


Substitution  of  (4.8)  for  the  lower  limit  in  (c.3)  results  in  its 
vanishing.  Hence,  considering  the  right  side  of  (C.3)  at  its  upper 
limit  as  t  v  results  in 


a 


V  Q  dr  =  lim 

vxv 

r-x» 

t-*v 


(rVt)(rQv)'  -  (rQv)(rVt)' 
v(v  +  1)  -  t(t  +  1) 


(C.4) 


The  numerator  vanishes  since  it  reduces  to  the  Wronskian  W  as  t  -*■  v. 

v 

This  also  implies  that 


V 

v 


a  Q 
v  v 


(C.5) 
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as  was  pointed  out  previously.  Then  applying  L*  Hospital’s  rule  to 
(C.4)  results  in 


/ 


dW 
_ t 

dt 

V  Q  dr  =  lim  .■  V*.. 

v  v  r^oo  +  !> 


(C.6) 


Since  the  Wronskain  of  any  two  independent  functions  is  a  constant  (C.6) 
reduces  to 

dW 


/ 


dt 


vvVr  "  -  T^TT) 


Making  use  of  (C.5),  we  have 


dW 
_ t 

dt 


t=v 


Vr) 

QV(D 


(2v  +  1) 


/ 


(Qv)  dr 


(C.7) 


And  substituting  (C.7)  in  the  expression  for  U(r,0),  we  have 


U(r 


Qv  0>)  Q^(t)  (cos  (tt  -  0)) 

♦0)  — 7  2^ 


(Qv(r))  dr 


sin  vtt 


for  r  _>  a. 

(C.8) 


The  integral  in  the  denominator  is  numerically  inconvenient .  To 
simplify  (C.8)  consider  the  expression  (C.3)  again 


VO  dr  = 
tv 


(rV  57  (rV  -  (rV  j?  (rV 

v(v  +  1)  -  t(t  +  1) 


(C.9) 


and  using  (C.5)  as 
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V 

t 


Mt 


then  the  normalized  integral  is 


dr  =  lim 
t+v 


Qt(r) 

Vfc(r) 


(rV  )  (rQ  )  -  (rQ  )  (r  Vj 

t  dr  v  xv  dr  t 

v(v  +  1)  -  t(t  +  1) 

(C.10) 


a 


Differentiating  numerator  and  denominator  with  respect  to  t,  and 
remembering  that  when  t  =  v,  an  eigenvalue,  the  function  V^(r)  is 
identical  with  Q^Cr). 

Hence , 


dr  -  lim 
t+v 


iF  <rVl 

If  (rV  -  (rQv)  If  <rVt>l 

a 

at 

v(v  +  1)  -  t(t  +  1)| 

(C.ll) 


The  denominator  of  (C.ll)  is  simply  -(2t  4-  1).  The  numerator  is 


lim 

t-n> 


<rV  iF  If  (rV  +  If  w»>  it  (tV 


<rV  it  If  <rV 


-  If  <rV  il  (rV 


(c.12) 


evaluated  at  r  =  a  and  r  =  00 .  Now  since  the  right  hand  side  of  (4.11) 
is  independent  of  t,  we  have 


__3 

at 


Sr  <*v 

(rVt) 


=  0 


r=a 


(C. 13) 


88 


or 


(rV  ft  dr  <rV 


'  W  <rV  37  (rV 


r=a 


=  0 


r=a 


(C. 14) 


or  since  when  t 


v  the  function  V  (r)  is  identical  with  Q  (r),  then 

V  V  5 


<rV  77  37  <rV 


-  37  (rV  37  W)  ‘  0 


(C.15) 


r=a 


Using  this  in  (C.12)  we  obtain  for  the  right  hand  side  of  (C.ll) 


/ 


(Q.)  dr  = 


(rV  ft  37  (rV  ~  37  (r<V  37  <rV 

~(2v  +  1) 


(C. 16) 


The  numerator  vanishes  as  t-**>  .  If  Q^Ca)  f  0,  that  is,  the  earth  is 
not  a  perfect  conductor,  then  we  can  write 


/ 


Q  (r)  dr  = 


-1 


(2v  +  1) 


<rv2  If 


Sf  (rQt} 

rQf 


r=a 

t=v 


(C. 17) 


If  the  earth  is  a  perfect  conductor  then  from  (4.10) 


(rVv) 


»  0 


r=a 


(C. 18) 


and  remembering  that  as  t-»-v,  we  have  from  (C.12)  that 
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-1 

(2v  +  1) 


3?  <rV  ST 


r=a 

t=v 


for  the  perfectly  conducting  earth. 


Using  (C. 17) 
dition  (4.11)  and 


and  (C.19)  in  the  expression  for  the  boundary  con 
(4.10),  we  obtain  expressions  (4.12)  and  (4.13). 


(C.19) 
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APPENDIX  D 


THE  MODIFIED  HANKEL  FUNCTIONS  AND  THEIR  PROPERTIES 


The  differential  equation 


d2u 

dz2 


+  z  u  =  0 


(D-l) 


is  known  as  Stokes1  equation  and  general  solutions  can  be  written  in 
terms  of  Bessel  functions  of  order  one-third  [22].  The  only  singularity 
of  this  equation  is  an  irregular  singularity  at  infinity.  Equation  (D-l) 
has  solutions  h^(z)  and  h^(z)  which  can  be  written  in  terms  of  the  two 
Bessel  functions  of  the  third  kind  of  order  one-third  as  follows: 


h^z)  =  (|  z3/2)1/3  (|  Z3/2) 


(D-2) 


h2(z)  =  (f  z3/2)1/3  H«>  (f  z3'2) 


(D-3) 


^1/3  anc*  ^1/3  are  tr^P-*-e  valued  functions  of  x,  with  a  branch 

point  at  the  origin.  However,  the  product  of  the  two  multipie-valued 
functions  (-j  z3^2)3^3  and  (-j  z3^2)  is  the  single  valued  function 

h . (z) .  The  Wronskian  of  h^  and  h2>  W(h^,  h^)  has  an  identically  van¬ 
ishing  derivative,  so  that 

4i  3 

w(hr  h2)  -  -  ~  (f) 

=  -1.  457  495  441  040  461  i  (D-4) 


a  constant. 
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The  solution  h^(z)  and  h^(z)  are  also  related  to  the  Airy  integrals 
Ai(z)  and  Bi(z),  defined  for  real  values  of  z  by 


Ai(-z)  =  —  /  cos  - zt)  dt 


(D-5) 


Then 


}i(_z)  =  y  [e_t 


-t3/3  -zt  .  .  ,t3  .  .  1 

+  sin  (-^ - zt)  dt 


(D-6) 


h1(z) 

Al(-z) 

-  i  Bi(-z)J 

(D-7) 

h2(z) 

■k‘! 

^Ai(-z) 

+  i  Bi(-z)J 

(0-8) 

Ai  (z) 

i 

2k 

h1(-z) 

+  *  h  (-z) 

(0-9) 

2k 

Bi  (z) 

i 

2k 

hj^C-z) 

-  h  (-z) 

2k 

(D-10) 

where 


or 


and 


o  2/3  r 

k  =  (|)  (i  -  H~) 


k  =  1.310  370  697  104  448  -  0.756  542  874  711  451  i 


2k 


=  0.286  178  560  638  333  +  0.165  225  269  020  841  i 


and  the  star  denotes  complex  conjugate. 
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Asymptotic  Expansions 


The  asymptotic  expansion  of  h^(z)  for  -  -j—  <  arg  z  <  , 


v,  t  ^  -1/A 

h^(z)~  a  z  e 


2/ 3i 


3/2  5ui  r 
12 


1  +  Y  (-!)“  C_  z“3m/2 


m=l 


m 


(9-A) (81-A)  ...  (9  [2m-l] 2  -A) 

.Am  _m 
2  3  m! 


(D-ll) 


— This  ray 
The  numerical 


a  =  21/3  31/6  tt  1/2  =  0.853  667  218  838  951. 


This  expression  cannot  be  used  along  the  ray,  arg  z  = 
is  the  branch-cut  for  the  multiple  valued  functions, 
coefficient  is 


An  asymptotic  expansion  valid  on  the  ray,  arg  z  =  -  ——  is 


o/o,  3/2  5rri 

-1/A  2/31  Z  12 

h1(z)~a  z  e 


[i  <- 

L  m=l 


-i)m  C  z_3m/2l 

m  J 


(D-12) 


-1/A 
+  a  z  e 


-2/3i  z 


3/2  _  llTTi 
12 


[l  +S 

L  m=l 


(i)m  C  z  3m/2 
m 


This  expansion  holds  for  -  —  <  arg  z  <  0;  the  branch  cut  for  the 

fractional  powers  of  z  can  be  drawn  anywhere  within  the  sector 

n  27T 

0  <  arg  z  <  — . 

The  existence  of  two  expressions  of  different  forms  which  represent 
asymptotically  the  same  integral  function,  h^(z),  is  an  example  of  Stokes’ 
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phenomenon . 


The  region  —  <  arg  z  <  —  for  (D-ll)  coincides  with  the  region 
4tt  J  2tt  J 

“  <  arg  z  <  -  -r—  for  (D-12).  Thus  in  this  region,  we  use  arg  z  =  a 

•J  -J  o  / 

in  (D-ll)  and  arg  z  =  a  -  2ir  in  (D-12)  with  y—  <  a  <  y— . 

The  asymptotic  expansions  for  h2(z)  are 


h2(z) 


a  z 


-1/4 


-2/3i  z 


3/2  5iri 
12 


m=l 


(i)m  C 


-3m/ 2 


m 


4tt 

valid  for  -  —  <  arg 


2tt 

z  <  — ,  and 


(D-13) 


9/-  3/2  5iTi 

.  ,  ,  -1/4  2/31  z  +  IT 

h2(z)  ~  a  z  e 


L  m=l 


(i)m  C  z'3m/2 
m 


(D-14) 


3/2  llTri 

-1/4  Z/J1  z  12  .  .  V'  ,  ,.m  _  -3m/2 

+  a  z  e  1  +  /i  (-i)  C  z 

i  m 

m=l 


valid  for  0  <  arg  z  <  — . 

Special  Properties  of  h^(z)  and  t^Cz) 


The  zeros  of  h1 (z)  and  those  of  h- 1 (z)  lie  on  the  ray  of  slope 

2tt  1  1 

-  —  ,  while  those  of  hjCz)  and  of  (z)  lie  on  the  ray  of  slope 

2tt 

angle  —  .  The  zeros  of  all  these  functions  are  simple. 

If 


T 


3/2 


a  /3  -n2/3 

and  z  =  (y  T) 


(D-15) 
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with 


Then 


and 


2 

arg  z  =  —  arg  T 

(T)  =  T_1/3  hx(z) 

(T)  =  T_1/3  h2(z) 

"2/3  m-<|)1/3  hl'U) 

1/2  +  ZjlL 

H2/3  <T>  “  <f>  e  3  T-2/3h2'(z) 


(D-16) 

(D-17) 

(D-18) 

(D-19) 

(D-20) 


Figures  D-l  and  D-2  give  the  regions  of  validity  of  the  asymptotic 
expansions  of  h^(z)  and  h„(z). 
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REGION  OF  VALIDITY 
OF  EQU.  (D-ll) 


REGION  OF  VALIDITY 
OF  EQU.  (D-12) 


Figure  D-l  .  REGIONS  OF  VALIDITY  OF  THE  ASYMPTOTIC  EXPANSION  FOR  h2  ( z ) 


REGION  OF  VALIDITY 
OF  EQU  (D— 13) 


REGION  OF  VALIDITY 
OF  EQU.  (D-l 4) 


Figure  D-2.  REGIONS  OF  VALIDITY  OF  THE  ASYMPTOTIC  EXPANSIONS  OF  h2(z) 
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